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We study a Kondo impurity model with additional uniaxial anisotropy D in a non-zero mag- 
netic field B using the Numerical Renormalization Group (NRG). The ratio g e /gs of electron and 
impurity g-factor is regarded as a free parameter and, in particular, the special cases of a "local" 
(g e = 0) and "bulk" (g e — gs) field are considered. For a bulk field, the relationship between the 
impurity magnetization M and the impurity contribution to the magnetization Mi mp is investigated. 
Furthermore, we study how the value of g e affects the impurity magnetization curves. In case of 
an isotropic impurity with g e = gs, it is demonstrated that at zero temperature Mi(B), unlike 
Mimp(B), does not display universal behavior. With additional "easy axis" anisotropy, the impurity 
magnetization is well described by a shifted and rescaled Brillouin function on energy scales that 
are small compared to \D\. In case of "hard axis" anisotropy, the magnetization curves can feature 
steps which are due to field-induced pseudo-spin- 1/2 Kondo effects. For large anisotropy and a 
local field, these screening effects are described by an anisotropic spin-1/2 Kondo model with an 
additional scattering term that is spin-dependent (in contrast to ordinary potential scattering). Our 
study is motivated by the question how the magnetic properties of a deposited magnetic molecule 
are modified by the interaction with a non-magnetic metallic surface. 

PACS numbers: 73.20.Hb, 75.30.Cr, 75.30.Gw, 75.50.Xx, 75.60. Ej 



I. INTRODUCTION 

Magnetic molecules offer the prospect of encoding and 
storing information in their magnetic state. The lat- 
ter point applies, in particular, to bistable molecules 
such as single molecule magnets (SMMs). The possi- 
bility to store, e.g., one bit of information in the state 
of a single molecule would constitute an enormous de- 
gree of miniaturization and could lead to data storage 
technologies with significantly increased areal densitypJ 
However, to make a (potentially elusive) technological 
application feasible, the molecules need to be individually 
addressable so that their magnetic state can be probed 
and manipulated on a molecule-by-molecule basis. In 
the last years, there has been an increasing interest in 
the question whether this functionality can be achieved 
by a controlled deposition of magnetic molecules on suit- 
able substratesP"^ While such an approach might solve 
the problem of addressability, it can introduce new com- 
plications due to interactions between the molecules and 
the surface. Depending on details such as the molecule's 
ligands, the presence of an additional decoupling layer, 
and, of course, the characteristics of the surface, the in- 
teraction with the substrate might alter the magnetic 
properties of the molecule in an important (and possi- 
bly adverse) way. Thus, even if the magnetic response 
of the isolated molecule is well understood (e.g., through 
a description by a suitable spin model-), its magnetic 
properties in contact with the surface have to be reinves- 
tigated. 

In this article, we study a single-channel Kondo impu- 
rity model with non-zero magnetic field and additional 
uniaxial anisotropy D(S Z ) 2 for the impurity spin oper- 
ator S. Such an anisotropy term (along with trans- 



verse anisotropy E[(S X ) 2 — (S y ) 2 ]) is a common part of a 
pure spin model for the description of isolated magnetic 
molecules (in particular, for representing SMMs) I ^ The 
quantum impurity model is intended to serve as a min- 
imal representation of an anisotropic magnetic molecule 
on a non-magnetic metallic substrate and, with trans- 
verse anisotropy E, has already been used to describe 
SMMs interacting with metallic electrodesP® Further- 
more, it has been found that the above uniaxial and 
transverse anisotropy terms are also appropriate to model 
the surface-induced anisotropy of a single magnetic atom 
on a metallic substrate with a decoupling layer To 
investigate how the interaction with the electrons affects 
the magnetic properties of the impurity, we carry out 
Numerical Renormalization GroupP^Hlll (NRG) calcula- 
tions and focus on the magnetic field dependence of the 
impurity magnetization. 

Regarding the experimental situation, the magnetic 
moment of deposited molecules (or atoms 'f^ can be mea- 
sured using metho ds s uch as X-ray magnetic circular 
dichroism (XMCD)P™ 21 XMCD is an element-specific 
technique of high sensitivity based on the absorption of 
circularly polarized X-rays and can be used to obtain 
an ensemble-averaged result for the magnetic field de- 
pendent molecule magnetization. In principle, it is also 
possible to extract information about different contribu- 
tions to the observed magnetic moment (such as the or- 
bital and spin contributio n) from the XMCD data using, 
e.g., sum rulesP ^ 18 * 20 * 23 lln the last years, magnetization 
curves of magnetic atoms on non-magnetic metallic sur- 
faces could also be recorded using spin-p olarized scanning 
tunneling spectroscopy (SP-STS)! 11 * 24 " 28 ! In contrast to 
XMCD, this method provides a time-average of the field- 
dependent magnetic moment of a single atom. It has 
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been demonstrated that SP-STS can also be applied to 
(suitable) deposited magnetic molecules! 29 " 31 ^ 

The static magnetization of Kondo impurity models 
(including related models such as the single- impurity An- 
derson model) has been investigated by a number of tech- 
niques. Among these are Green's-function methods ] 32 * 33 ^ 
the Bethe AnsatzJ^Sl an d NRQ^HH (including density 
matrix based extensions) . By now, there are also several 
studies of the time dependence of the magnetization in 
non-equilibrium situations (e.g., after a quantum quench 
or with a non-zero voltage bias)P3EH In particular, non- 
equilibrium spin dynamics of impurity models can be in- 
vestigated by using a genera lization of NRG called time- 
dependent NRG (TD-NRG) PSM] 

The present article extends existing NRG results for 
the Kondo model with uniaxial anisotropy^ to the case 
of non-zero magnetic field. The system with non-zero 
field (with a focus on the properties of spectral functions) 
has been previously studied in Refs. [571 and [551 Further- 
more, magnetization curves for isotropic and anisotropic 
Kondo impurities have been calculated in Ref. HgJ We 
would like to stress, however, that our investigation 
places emphasis on different aspects of the problem and 
is thus complementary to Ref. HU 

The remainder of this article is organized as follows. 
In Sec. |ITJ the quantum impurity model is introduced 
and transformed to a representation that is suitable for 
further numerical treatment. Sec. Ill provides infor- 



mation about our use of the NRG method and contains 
definitions of the considered observables. In Sec. |IV| we 
study the relationship between the impurity magnetiza- 
tion and the impurity contribution to the magnetization 
for an isotropic system (D = 0) . After an investigation of 
the Kondo model with additional "easy axis" anisotropy 
(D < 0) in Sec. |vj the case of "hard axis" anisotropy 
(D > 0) is considered in Sec. |VI| There, we also clar- 
ify how a non-zero magnetic coupling of the conduction 
electrons affects the impurity magnetization. Further- 
more, an effective model is derived in order to describe 
the field-induced pseudo-spin- 1/2 Kondo effects that are 
observed in the magnetization curves for large D. We 
conclude this article with a summary of the results in 



Sec. VII and a brief description of the technical details 
of an NRG calculation with non-zero magnetic coupling 
of the conduction electrons in App. [AJ 



II. MODEL 



A. Hamiltonian 



In this work, we study a Hamilton operator if consist- 
ing of three parts: 



H — if electrons + H 



coupling ~t~ if impurity 



(1) 



The first term, if electrons, represents non-interacting 
tight-binding electrons whose hopping between different 



sites i and j of a periodic lattice with L sites is described 
by the corresponding hopping parameter tif 



H electIons = 53 t 



rft d 



(2) 



Here, d^J is a destruction (creation) operator for an elec- 
tron with spin projection a = ±1/2 = j" / \, at lat- 
tice site i. The effect of an external magnetic field B 
is taken into account by a Zeeman term with electron 
g- factor g e , Bohr magneton /ib, and the z-componcnt 
of the total spin of the electrons S z = \ Si fet — 24) 

with n icr = d ia di a . Using a discrete Fourier transforma- 
tion, c' ka = (1/vT) J2j e%k Hj £^cr' Hamiltonian (2| can 
be cquivalently written in the more common form 



h.a 



£fc + crgeHBB) cl a c k<T 



(3) 



= e k „(B) 



with a dispersion relation Eka(B), assigning an energy 
e to a wavevector fc, that now depends on spin projec- 
tion and magnetic field. In general, the spin-independent 
dispersion relation e k is anisotropic in fe-space. 

For the interaction term in Eq. ([I]) , we use a standard 
isotropic Kondo coupling, 



H 



coupling 



= JS ■ s , 



(4) 



and assume that the impurity spin S couples antiferro- 
magnetically ( J > 0) to the electronic spin at the origin, 
which is given by s = (1/2L) Y,k,k>,w £fc / i V£fe'" with 
the vector of Pauli matrices <j . 

Finally, the impurity part of Hamiltonian ([I]) repre- 
sents a localized spin with quantum number S which 
couples to the external magnetic field with g-factor gs 
and possesses an additional uniaxial anisotropy D: 



-ff impurity — D(S Z ) 2 + gsLLsBS* 



(5) 



With the chosen convention, the impurity spin has 
an "easy axis" for D < and a "hard axis" or 
an "easy plane" for D > 0. A further transverse 
anisotropy E[(S X ) 2 — (S v ) 2 ] is not considered in this ar- 
ticle, if impurity can be seen as a minimal representation 
of a magnetic molecule with a single magnet ic ce nter or 
as a "giant spin approximation" for an SMMP^H 

Hamiltonian ([T]) corresponds to an anisotropic single- 
channel Kondo impurity model in an external magnetic 
field. The special choices g e = and g e — gs for the 
electron g-factor are referred to as a "local" and "bulk" 
magnetic field, respectively. Regarding the modeling of 
a deposited magnetic molecule, it has to be emphasized 
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that the Hamiltonian suffers from a number of simplifica- 
tions. For example, there is no orbital contribution to the 
magnetism, and no charge fluctuations between molecule 
and surface are possible. In this article, we only con- 
sider the effect of the Kondo coupling on the magnetic 
properties of the impurity spin. 



B. Transformation to an energy representation 

In order to treat Hamiltonian |TJ) using NRG, i? e iectrons 
and -ff coupling are expressed via a continuous energy rep- 
resentation for the electronic degrees of freedom. To 
this end, we first take a standard continuum limit in 
fc-space OLe., we consider a lattice of dimension d with 
L ^> l)p31 By adapting the corresponding expression 
for the two-impurity Kondo model from Ref. [60] to the 
single- impurity case (see also Ref. loT]) . we then define 
those states with energy e to which the localized spin 
directly couples: 



V(2tt)<V(£ - W) 



/"dfe(5(e-e M (fe,S)) c feM , (6) 



where we have introduced the abbreviation h — g e pBB 
and the normalized density of states (DOS) per spin pro- 
jection and lattice site p(e) = (l/L)J2 k S(e — De- 
noting the half- width of the conduction band by W, the 
allowed energies e for spin projection p span the interval 
[— W + ph, W + ph]. The new operator a efl is properly 
normalized because of the pre-factor involving the DOS. 

If we are only interested in impurity properties, then 
all other electronic states different from those defined in 
Eq. © can be safely dropped without introducing any 
approximation.^ This leads to the desired continuous 
energy representation of Hamiltonian : 



-W+ph 

+ J § J2\ &e^p(e-ph) at p 



/W+vh \ 
de' y/p{e' - vh) a e , v + H 
-W+vh I 



impurity 



(7) 



For h = 0, i.e., for B = or g e = 0, Eq. ([7| reduces to 
the well-known expression for the energy representation 
of the Kondo modelP^l In the following, we consider the 
case of a constant DOS: p(e) — 1/2W — p. 



III. METHOD AND OBSERVABLES 
A. Method: NRG 

Approximate eigenvalues and eigenvectors of Hamil- 
tonian |7]) for the calculation of impurity properties 
can be obtained using the Numerical Renormalization 
Q r0U p(TlHlll (NRG). However, the procedure leading to 
the parameters of the Wilson chain has to be slightly 
modified if h ^ (see App. [A] for a brief discussion of 
the required changes). 

A non-zero magnetic field breaks the full SU(2)- 
symmetry in spin space of Hamiltonian ([T]). For this 
reason, we label eigenstates of H only with the charge 
quantum number Q and the magnetic quantum number 
Stotai of the z-component of the total spin. Except for 
one example in Sec. |VI[ all NRG calculations are carried 
out using th e dis cretization scheme proposed by Zitko 
and PruschkePEU with averaging over 4 z-values that are 
equidistantly spaced on the interval (0,1]. Observables 
are computed using only states that are kept after trun- 
cation and results are averaged over even and odd sites of 
the Wilson chain according to the prescription of Ref. HU 
We use a discretization parameter A = 3, a dimensionless 
inverse temperature /? = 0.7, and a fixed number of kept 
states of the order of N = 5000 to achieve convergence 
for all considered observables within the resolution of the 
presented plots. Nevertheless, at A > 1 there might still 
be slight systematic deviations for non-zero temperature, 
which can for example be demonstrated by setting J = 
and comparing the NRG results with the analytical solu- 
tion for a free spin. It is necessary to perform a separate 
NRG calculation for each value of the magnetic field. If 
curves are shown in a plot, they are thus the result of 
a spline interpolation through the numerically obtained 
data points. 



B. Observables 

In our calculations we focus on the impurity magneti- 
zation which is defined as the expectation value of the im- 
purity magnetization operator with respect to the eigen- 
states of the total Hamiltonian (JlJ: 



M(T, B) = - 



^-^impurity 

dB 



-9sHb{S z 



'total 



total 



(8) 

Furthermore, we consider the impurity contribution to 
the entropy, magnetization, and magnetic susceptibility. 
The impurity contribution to some quantity O is defined 
in the usual wayP^ 



o 



with impurity 
total 



o 



w/o impurity 
total 



(9) 
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The observable O 



w/o impurity 
total 



for the system without im- 



purity is also calculated using NRG by removing the 
impurity part from the Wilson chain. For the en- 
tropy S(T,B), the magnetization M(T,B), and the 
susceptibility x(T,B), we use the standard definitions 
S(T,B) = -dQ(T,B)/dT, M(T, B) = -dtt(T,B)/8B, 
and x(T,B) = dM(T, B)/dB, with n(T,B) being the 
grand-canonical potential. If the electron g-factor is zero, 
we have the special case M(T, B) = Mi mp (T, B). 

In the grand-canonical calculations the chemical po- 
tential is assumed to be zero. For a symmetric DOS, 
p(e) — p{— e), the free electron band is thus on average 
half-filled for arbitrary magnetic field and temperature. 



IV. ISOTROPIC IMPURITIES 

Let us first consider the case of an isotropic impurity 
with D = in Hamiltonian ([5| and study the impurity 
contribution to the magnetization M- lmv and the impu- 
rity magnetization Ai, both as function of temperature 
and magnetic field. For the moment, we are only con- 
cerned with the special case of equal g-factors of impu- 
rity and electrons (corresponding to the case of a bulk 
magnetic field). Recalling the motivation given in the in- 
troduction, Ai as the expectation value of the impurity 
magnetization operator should be the observable that is 
more closely related to experimental magnetization data 
obtained by methods such as XMCD. Note that M. and 
A/imp become equivalent if impurity and electrons decou- 
ple (which happens for J — > or T — > oo). 



A. Field dependence of the magnetization 

In case of the isotropic Kondo model with g e = 
gs and arbitrary impurity spin S, the Bethe Ansatz 
(BA) allows for the derivation of a closed expression for 
the impurity c ontribut ion to the magnetization at zero 
temperature.^" ^ 7 * 39 * 40 * M; mp is known to display univer- 
sal behavior in the so-called scaling regime, in which all 
relevant energy scales are small compared to the energy 
cutoff (or the finite bandwidth) Bare parameters of the 
model can then be absorbed into an energy scale k^Tfj so 
that the field-dependence of M- lmp at T = is described 
by a universal function f(x), with x being the rescaled 
magnetic field: x = gspsB /IcbTh^Th is related to the 
Kondo temperature Tjc by a constant factorP^ 



H 



K 



(10) 



In Fig. [T]wc plot the Bethe Ansatz solution for M; mp 
from Ref. [591 for three different impurity spins S = 1/2, 1, 
and 3/2.^ J(x) is a strictly monotonically increasing 
function of x and approaches the saturation magneti- 
zation of a free spin, gsPB>5, f or x — > oo with slowly 
decaying logarithmic corrections! 39 * 40 * 6 -^ The behavior in 
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FIG. 1. (Color online) Main plots: Impurity contribution 
to the magnetization il/i mp and impurity magnetization M 
as function of magnetic field for g e = gs, three different cou- 
plings p,J, and for impurity spin a) S = 1/2, b) S = 1, and c) 
S = 3/2. The temperature is k B T/W w 1.54 ■ 10" 15 w and 
the field is rescaled using ksTn- For the shown A4 -curves 
the coupling pj always increases from bottom to top. The 
universal BA solution at T = has only been calculated for 
gspBB > ksTn and gspsB <C fesTf/ in plots b) and c). Up- 
per left insets show NRG results for Mi mp at T ~ (black 
lines) and finite temperature (solid gray lines) as a function of 
magnetic field, now expressed in units of W. M- lrnp for J = 
is also computed using NRG and resembles the magnetization 
of the free spin. Thermal energies increase from left to right 
and range from 1.79- 10~ 6 W (plot a)) or 1.95 ■ 10" 12 1F (plots 
b) and c)) to 6.79- 10 _3 VK. Lower right insets show a close- 
up of the magnetization curves for b) low fields and c) high 
fields along with data points for Mj mp that are multiplied by 
a constant > 1. 
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the limit x —> depends on the value of 5: In case 
of 5 = 1/2, f(x) oc x for small x, whereas for 5 > 1 
the function f(x) goes to the saturation magnetization 
of a reduced s pin with 5 — 1/2, again with logarith- 
mic corrections! 39 * 40 * 6 ^* This low-field behavior mirrors 
the Kondo screening which, for vanishing magnetic field, 
reduces the impurit y spin 5 to a residual spin S — 1/2 in 
the limit T/T K < lP^The magnetic properties of the 
impurity are furthermore markedly different from that of 
a free spin as the magnetization of a free spin at T = 
saturates for any non-zero magnetic field. 

Using NRG, we have calculated Mj mp for several values 
of the coupling strength p J and have fitted the obtained 
curves to the respective universal BA curve by employing 
Tfj as a fit parameter (see Fig. [I]). The nice agreement 
with the Bethe Ansatz solution demonstrates the univer- 
sal behavior that M- lmv displays for small pJ and allows 
us to reliably determine the value of T B even for impurity 
spin S > 1. Note that for very large magnetic fields, i.e., 
for gsp B B < W, we leave the scaling regime and Mi mp , 
as calculated by NRG, starts to drop below the universal 
BA curve (this is not shown in Fig. [TJ. The determined 
approximate values of k B T B /W are found in Table [Ij In 
the cutoff scheme used in the BA solution of the Kondo 
model with arbitrary impurity spin, the dynamically gen- 
erated low-energy scale T , whose ratio to T B is a uni- 
versal number, d epends on the value of 5 (before taking 
the scaling limit) F^ mD Herej we find th&t thc fittcd 

values of Th increase with the impurity spin for fixed 
coupling strength and, furthermore, that the relative de- 
viation between the results for different 5 decreases when 
pJ is reduced. However, even for the smallest consid- 
ered coupling strength (pJ — 0.05), the values of Th 
for 5=1/2 and 5 = 3/2 still deviate by about 44 %. 
According to Eq. (fIo| ), the values of T H for 5 = 1/2 
reported in Table [if correspond to the following Kondo 
temperatures: k B T K /W « 4.79 • 10" 10 {pJ = 0.05), 
1.80 • 10~ 7 (0.07), and 5.08 • 10~ 6 (0.09). These re sults 
can be compared with the standard estimate for Tr-* 12 * 13 * 

k B T K /W « v/^Jexp (-1/pJ) , (11) 

which is valid for small coupling and gives k B Tx/W ps 
4.61 • 10~ 10 (pJ = 0.05), 1.65 • 10" 7 (0.07), and 4.48 • 
10~ 6 (0.09). As a further check, we have determined the 
Kondo temperature for 5 = 1/2 and pJ = 0.07 by fit- 
ting the BA solution for the impurity contribution to the 
susceptibility from Ref. [39] and the impurity contribu- 
tion to the entropy from Ref. This gives a value of 
k B Tx/W « 1.79 • 10~ 7 , which is quite similar to the one 
following from Table [I] 

The upper left insets of Fig. [T] show finite temper- 
ature NRG results for Af; mp with a coupling strength 
pJ = 0.07. While the Bethe Ansatz provides a closed 
expression for Mi mp at zero temperature, a calculation 
for finite temperature leads to so-called thermodynamic 
BA equatio ns th at, at least in general, have to be solved 
numerically! 38 * 70 ^ Hence, finite temperature results for 



the magnetization are not easily available. As a refer- 
ence point, we replot the zero temperature magnetization 
curve that crosses over to the strong coupling regime in 
the vicinity of gsp^ B B ps k B T B - As long as the ther- 
mal energy is small compared to the Zeeman energy, the 
magnetization always closely follows the zero tempera- 
ture curve. On the other hand, if the thermal energy is 
not negligibly small compared to the Zeeman energy, we 
have to distinguish between complete screening and un- 
derscreening of the impurity spin. For 5 > 1, finite tem- 
perature is always important as it also affects the residual 
spin. On the energy scale gsp B B ~ k B T there is a swift 
drop of M 1Ta p that is eventually followed by a linear de- 
cay for small fields gs[i B B <C k B T. In the special case 
5 = 1/2, however, finite temperature has little effect if 
T <C Tk and the magnetization already displays a linear 
dependence on the magnetic field for gsp B B ps k B T due 
to the Kondo screening. In the upper left insets of Fig. 
[T] we also compare the results for M imp with NRG cal- 
culations for vanishing coupling J — 0. This comparison 
is meant to illustrate the influence of a non-zero value of 
Jp3 At high temperatures (compared to the Kondo tem- 
perature) the impurity is progressively decoupled from 
the electronic system and its magnetization will hence 
resemble the result for J = more closely. However, 
note that the impurity only becomes asymptotically free 
for high temperatures. 

In addition to the impurity contribution to the magne- 
tization -Mi mp , we also plot the impurity magnetization 
M in Fig. [I] for the same values of the coupling pJ 
and negligible temperature. The magnetic field is again 
rescaled by k B T B using the values from Table [TJ We find 
that M. and Mi mp differ for all considered magnetic fields 
with M. being larger than Mj mp . This means, in partic- 
ular, that for large magnetic fields M. comes closer to 
the saturation magnetization of a free spin than M lmp 
does and, furthermore, that the magnetization of the con- 
duction electrons is reduced due to the interaction with 
the impurity spin. Upon decreasing pJ we observe that 
the impurity magnetization becomes smaller and thus ap- 
proaches the universal curve for Mi mp . A comparison of 
the NRG results for M. and M imp shows that both quan- 
tities are proportional to each other, i.e., A4 — aAl imp 
with a proportionality factor a > 1 that depends on the 
coupling strength p J (see Table [I] for a list of the com- 
puted values of a). This proportionality is illustrated for 
the case of small magnetic fields (for 5=1) and large 
magnetic fields (for 5 = 3/2) in the lower right insets of 
Fig. [l] While the obtained values for a decrease with 
increasing impurity spin 5, the values for different 5 dif- 
fer by less than 0.2 % according to Table [TJ Since impu- 
rity and electrons progressively decouple at high temper- 
atures, we expect a to be temperature dependent with 
a — >• 1 for k B T/W ^> 1. The results presented in Fig. [I] 
show that the magnetic field cannot be rescaled by k B T B 
or any scale proportional to it so as to produce a universal 
curve for the impurity magnetization Ai. 

To elucidate our findings we refer to one of the orig- 
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TABLE I. Approximate values of UbTh /W as used in Fig. [Tj obtained by fitting the universal Bethe Ansatz solution for Mi mp , 
and proportionality factors a(pj) relating M and Mi mp according to M — a?Mj mp . The results for a have been averaged over 
magnetic fields qs^bB /W G [10 -13 , 10 _1 ] for ksT/W fa 1.54 ■ 10 - ~ 0. Numbers in parentheses give the corresponding 
standard deviation for the last decimal place. For Zeeman energies close to the band edge (i.e., gsHBB < W), which are not 
considered for the average, a noticeably decreases (increases) for S = 1/2 (S = 1,3/2). 





S = 


1/2 




S = 1 


S = 


3/2 


pj 


k B T H /W 


a 


k B T„/W 


a 


k B T H /W 


a 


0.05 


7.29 ■ W~ w 


1.02659(1) 


8.49 • 10" iu 


1.026503(7) 


1.05- 10" 9 


1.02638(2) 


0.07 


2.74- 10 -7 


1.03822(2) 


3.39- 10~ 7 


1.03792(3) 


4.55- 10~ 7 


1.03751(6) 


0.09 


7.72 ■ 10~ 6 


1.05048(3) 


1.02 • 10~ 5 


1.04970(8) 


1.51 • 10~ 5 


1.0487(2) 



inal Bethe Ansatz investigations of the Kondo modelP^ 
There it is found that, under the assumptions of a BA 
calculation (including an arbitrarily large energy cutoff 
@), M. = Mi ro p. To study the influence of the cutoff 
scheme, a comparison with perturbation theory is car- 
ried out showing that A4 has leading corrections of the 
order l/ln(£^), whereas the corrections of M lmp vanish 
like l/2> and thus much faster P2l The regime in which 
all relevant energy scales are negligible compared to the 
cutoff in a logarithmic sense, e.g. h\{S> / qs^bB) 1, is 
termed "extreme scaling limit" P2 

With this background we reach the following inter- 
pretation of our NRG results for M. and M; mp : For 
the chosen values of the coupling strength pj the half- 
bandwidth W (basically serving as the unit of energy) 
can be regarded as very large compared to all relevant 
energy scales £ so that corrections of order £ /W can be 
expected to be small. It is for this reason that we find 
nice agreement with the universal BA solution for M- lmp . 
On the other hand, corrections of order l/\n(W/£) are 
not necessarily negligible for a finite value of W. This ap- 
pears to be an adequate explanation for our NRG results 
showing that M ^ M imp . Moreover, a decrease of pj 
corresponds to an increase of the bandwidth and thus 
bandwidth-related corrections should become smaller. 
Accordingly, M. decreases for smaller coupling strength 
and approaches M; mp . These observations might also 
bear some importance for experimental situations: While 
experimental parameters are certainly suitable to con- 
sider the scaling regime (in case the system exhibits uni- 
versal behavior), it is less clear whether an experimental 
system can be placed in the extreme scaling regime. 
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FIG. 2. (Color online) Saturation value in the limit T — s> of 

a) the impurity contribution to the magnetization Mi mp and 

b) the impurity magnetization M as function of the coupling 
strength pj for impurity spin S = 1, g e ~ gs, and several 
magnetic field values. The vertical lines mark a coupling of 
pj' w 0.276 for which we have checked by comparing with the 
Bethe Ansatz solution shown in Fig. [TJthat Mj mp still exhibits 
universal behavior (a fit of the BA curve gives ksTn /W ~ 
8.85 • 10 -2 ). Remaining lines are intended as a guide to the 
eye. 



B. Dependence of the zero temperature 
magnetization on the coupling strength 

To further illustrate the difference between the impu- 
rity contribution to the magnetization and the impurity 
magnetization, we examine how both quantities depend 
on the coupling strength p J for non-zero magnetic field at 
zero temperature. As before, the case of equal g-factors 
for impurity and electrons is considered. NRG results for 
impurity spin S = 1 and S = 3/2 are shown in Fig. [2] 
and Fig. [3j respectively. 



Let us begin the interpretation of the plots by consid- 
ering the limiting cases pJ — > and pJ — > oo. For vanish- 
ing coupling strength J — ¥ 0, impurity and electrons are 
decoupled and thus both Mi mp and M correspond to the 
magnetization of a free spin which saturates for any non- 
zero value of the magnetic field at T — 0. However, the 
behavior in the limit J — > oo differs for the two quantities, 
demonstrating that in general Mj mp ^ M. The values 
of Mj mp and M in this limit can be understood by con- 
sidering a simplified model: In the interaction term Eq. 
(|4|), the spin operator s depends on the occupation of 
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FIG. 3. (Color online) Saturation value in the limit T — > of 
a) the impurity contribution to the magnetization and b) the 
impurity magnetization as function of the coupling strength 
pj as in Fig. [5J but here for impurity spin S = 3/2. Mi mp 
again shows universal behavior for pj' ~ 0.276 and a fit of 
the universal BA curve now gives ksTn /W ~ 3.05 • 10 _1 . 



the lattice site that it is associated with. For very large 
values of J it is energetically favorable that this lattice 
site is singly occupied with a probability near one. We 
can then replace So by a spin-1/2 operator s in Eq. (Ej). 
Furthermore, all other terms in Hamiltonian that in- 
volve degrees of freedom different from the impurity and 
the lattice site to which it couples can be neglected. We 
are then left with a strongly coupled antifcrromagnetic 
dimer in a magnetic field with Hamiltonian 



H 



J— >oo 



= JS-s+g s p: B B(S z + s z ) 



(12) 



To determine the limiting value of AA/gs^B for pj — > 
oo at any positive magnetic field, we thus have to calcu- 
late the expectation value of — S z with respect to that 
eigenstate of total spin which has the lowest value of 
Stotai an d corresponding z-projection M tota i = -5 to tai- 



Denoting eigenstates of the effective Hamiltonian ( 12 ) as 
I Stotai, -Mtotai), we find: 



(0,0|5 Z |0,0) =0 for 5=1/2, (13) 



1 1 

2'~2 



s r - 



1 1 

2'~2 



(1,-115*11,-1) 



for 5=1, (14) 
for 5 = 3/2. (15) 



In contrast, in the limit pj —> oo the impurity con- 
tribution to the magnetization Mi mp /gsPB reduces to 
— (S z + s z ) with respect to the above eigenstates of 
Hamiltonian (12). This expectation value gives 5— 1/2. 



The case 5 = 1/2 is therefore special since both M llnp 
and M. go to zero for pj — s- oo. Figs. [2] and [3] show that 
NRG, as a method that is non-perturbative in J, can 
in fact reproduce the limiting values for large coupling 
strength. 72 

Let us now consider the magnetization for interme- 
diate values of pj. In the special case 5 = 1/2 both 
Mimp and M. are monotonically decreasing functions of 
the coupling strength that show similar behavior. We 
furthermore find that a larger value of the magnetic field 
B also leads to a larger value of the magnetization. In 
fact, all magnetization curves that we have calculated in- 
crease monotonically with B. On the other hand, -Mi mp 
and M. display qualitatively different behavior for im- 
purity spin 5 > 1. While M- lmp is again a monotonically 
decreasing function of the coupling strength, M develops 
a minimum for all considered magnetic fields. 

It is instructive to compare the NRG results with the 
Bethe Ansatz solution for Mi mp . According to the stan- 



dard estimate for the Kondo temperature from Eq. (Ill, 
ksTx/W is a monotonically increasing function of the 
coupling strength pj . Furthermore, a larger magnetic 
field always results in a larger magnetization for the sys- 
tem under consideration. For this reason, an increase of 
the coupling p J leads to a larger Kondo temperature Tk 
and thus, after rescaling, to a lower effective field and 
hence a smaller value of M imp . The impurity contribu- 
tion to the magnetization should therefore be a monoton- 
ically decreasing function of the coupling strength in the 
scaling regime. This observation is consistent with our 
NRG results for Mi mp as reported in Figs. [2] and [3] The 
vertical lines in the plots mark a coupling strength pj' 
for which we have checked by comparing with the BA 
solution that Mj mp still displays universal behavior. It 
can thus be seen as a lower bound for the couplings that 
can be assigned to the scaling regime. In contrast, the 
impurity magnetization Ai has a minimum for coupling 
strengths that are smaller than pj' . We conclude that 
this behavior of M. is not compatible with the standard 
scaling picture as described above. 



V. IMPURITIES WITH EASY AXIS 
ANISOTROPY 

We now deal with the case of an impurity with addi- 
tional easy axis anisotropy (i.e., with anisotropy param- 
eter D < in Eq. ([5])). In this section, emphasis is 
placed on the field dependence of the impurity magneti- 
zation Ai, again for the case of equal g-factors. Before 
considering the full impurity model given by Eq. Q, let 
us briefly recapitulate the magnetic properties of a free 
anisotropic spin with Hamiltonian (j5j). 

For negative anisotropy parameter D and vanishing 



magnetic field, the groundstate of a spin 5 > 1 is a dou- 
blet composed of the states with magnetic quantum num- 
ber M — ±5. In the special case 5 = 1/2, the anisotropy 
term D(S Z ) 2 evaluates to a constant and is thus insignif- 
icant for the thermodynamics. The first excited state 
is a singlet with M = for 5 — 1 and a doublet with 
AI = ±(5 — 1) for all larger spins. It follows that the 
energy gap between groundstate and first excited state is 
given by \D\(2S — 1). For thermal energies that are small 
compared to this gap, the zero-field magnetic susceptibil- 
ity approximately obeys a Curie law with Curie constant 
((5 Z ) 2 ) = S 2 (instead of 5(5+l)/3 for an isotropic spin). 

What do we expect for the full impurity model if there 
is an additional easy axis anisotropy? Since the ground- 
state doublet of the free anisotropic spin has AM = 25 > 
1, the two states it is comprised of are not connected by 
a single spin- flip, which changes M by 1. Furthermore, 
for increasing values of \D\ the gap in the energy spec- 
trum of the free anisotropic spin progressively suppresses 
scattering processes connecting groundstate and first ex- 
cited state. With the scattering picture in mind, one 
would thus assume that the Kondo effect is weakened 
by a negative value of D. This is in line with the simpli- 
fied picture for the limit \D\ — > oo: The anisotropy term 
D(S Z ) 2 then effectively acts as a projection operator onto 
the groundstate doublet of the impurity with M = ±5 
and hence asymptotically reduces the full Kondo interac- 
tion of Eq. Q to an Ising-type coupling.^ With respect 
to the impurity magnetization A4, there appears to be 
an even simpler argument: A larger negative value of 
the anisotropy parameter D energetically lifts all excited 
states of the impurity, which have reduced magnetic mo- 
ment in comparison to the groundstate doublet. At large 
|£>| one would thus expect that the excited states have 
less weight in the many-body groundstate of the full im- 
purity model leading to an increased value of M at zero 
temperature for positive magnetic field. 



A. Field dependence of the impurity magnetization 
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FIG. 4. (Color online) Impurity magnetization M for differ- 
ent anisotropy parameters D < (easy axis anisotropy) as 
function of magnetic field for k B T/W w 1.54 • 10" 15 « 0, 
coupling strength pj = 0.07, and impurity spin a) S = 1, b) 
S = 3/2, and c) S = 2. For the solid curves the value of \D\ in- 
creases from bottom to top and we compare with the magneti- 
zation of an isotropic impurity (dashed line). As before, equal 
g-factors of electrons and impurity are assumed. Vertical lines 
mark the value of ksTn /W , which is determined by fitting 
the universal Bethe Ansatz solution for M lmp with D = 0. 
For P J = 0.07 and S = 2, we find k B T„ /W « 6.8 • 10 -7 . 



In Fig. [4] low-temperature NRG results for the impu- 
rity magnetization M(B) for impurity spin 5 = 1, 3/2, 2 
are presented. We start the discussion of the results at 
high magnetic fields and move from there to lower fields. 
If the Zeeman energy is much larger than the anisotropy 
parameter, i.e., if gsl^BB 3> \D\, nearly isotropic behav- 
ior of M is observed. At smaller fields gs^BB « \D\, the 
impurity magnetization for D < begins to deviate from 
the isotropic curve and, for gs(J.BB <C \D\, converges to 
a D-dcpcndent value larger than gs/is(5 — 1/2). In the 
limit of low fields, the impurity magnetization curves for 
D < shown in Fig. [4] are well described by a linear law: 

M{B) w M {a) (D) + 7(D) • gsHsB/W . (16) 

A4^ \D) thus corresponds to the impurity magnetization 
in the groundstate of Hamiltonian ([l]) for infinitesimal 



magnetic field. The low-field behavior for D < as de- 
scribed by Eq. (16 1 is different from that displayed by 
an isotropic impurity: For D = and 5 > 1, the im- 
purity contribution to the magnetization Mj mp , which is 
proportional to M. in the isotropic case as demonstrated 
in Sec. |IV A| approaches the limit of zero m agnetic field 
with slowly decaying logarithmic corrections! 39 * 40 * 6 -^ 

From the results presented in Fig. [4] we conclude that 
for non-zero magnetic field and D < a larger value of 
\D\ leads to a larger impurity magnetization Ai, with an 
upper bound for \D\ — > 00 given by the free saturation 
value gsfJ-BS. This observation is in agreement with the 
expectations formulated at the beginning of this section. 
One might therefore say that an easy axis anisotropy sta- 
bilizes the impurity spin. 

Taking another look at Fig. [4] and focusing on the 
regime of small magnetic fields with Zeeman energy 
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FIG. 5. (Color online) Impurity magnetization M as function 
of magnetic field for fixed anisotropy D/W = — 10~ 3 and 
temperature k B T/W w 1.03 • 10" 6 for 5 = 1, k B T/W w 
1.28- 10 -8 for S = 3/2, and k B T/W w 1.58 • 1GT 10 for S = 2. 
The impurity spin 5 increases from bottom to top. Note that 
Ml is not saturated for any field in the plot range (cf. Fig. 
[4|. Open symbols represent fits using a rescaled and shifted 
Brillouin function, f(x) = 'yBsijjx), and solid (green) lines 
fits using a rescaled and shifted Langevin function, g(x) = 
jL(r}x). 



gsHsB <C \D\, one might be misled to think that there 
is a saturation of the impurity magnetization M. (this 
impression would not occur for an isotropic impurity). 
This raises the question whether it is possible to approx- 
imately describe the field dependence of M for gs^BB <C 
\D\ using a model for a free spin. In the simplest case, 
such a description could be provided by a Brillouin func- 
tion Bg{x), which gives the temperature and field depen- 
dence of the magnetization of a free and isotropic spin S. 
As demonstrated in Fig. [3] for one value of D, it is in fact 
possible to adequately fit the NRG results for M from 
Fig. [4] using a rescaled and shifted Brillouin function, 
f(x) = ~/Bs(r)x) with free parameters 7 and ry, as long 
as ksT <C \D\ and gs^sB <C \D\. However, for larger 
fields the behavior of the impurity magnetization clearly 
differs from the prediction of a Brillouin function, as seen 
in Fig. [4j The ratio of \D\ and ksTn (or alternatively 
^bTk) determines the "apparent saturation value" of M. 
and thus the parameter 7. In contrast to a fit with a 
modified Brillouin function, a classical description using 
a rescaled and shifted Langevin function, gix) = jL(rjx) 
with L(x) = lims_ i . 00 Bs{x), does not work well for mag- 
netic fields close to the "saturation field" (cf. the solid 
lines in Fig. [5]), as is to be expected for a quantum me- 
chanical system with low spin. Nevertheless, a fit using 
g(x) can produce reasonable results for fields that are 
small compared to the "saturation field" . 

The results depicted in Fig. [5] might be of importance 
for an experimental study of a system that is (approxi- 
mately) described by Hamiltonian with a strong easy 
axis anisotropy. It is then conceivable that a measure- 
ment of the magnetization for magnetic fields that can 



be realistically produced in an experiment (depending 
on the value of D, fields with gs^BB ss \D\ might not be 
obtainable) does not allow to distinguish between the be- 
havior of an anisotropic impurity spin and that of a free 
spin. Such a scenario seems more likely if the experimen- 
tal control over the g-factor and the absolute magnitude 
of the magnetization is limited, and if \D\ is large com- 
pared to ksTx so that the "apparent saturation value" 
of the impurity magnetization lies close to the free satu- 
ration value gsfJ-BS. 



B. 



Impurity contribution to the magnetization and 
the susceptibility 



As for an isotropic impurity (cf. Sec. IV A), we have 
analyzed the connection between the impurity magneti- 
zation M, which is shown in Fig. [4j and the impurity 
contribution to the magnetization Mj mp (not shown) for 
D < 0. It is found that the relation between both quanti- 
ties is the same as in the isotropic case, i.e., A4 = aM[ mp 
with a proportionality factor a that is independent of 
the anisotropy parameter D when taking into account 
the precision of the results from Table [TJ For impurity 
spin S = 2 and coupling strength pj = 0.07 the obtained 
value of the proportionality factor is a = 1.03701(3). 

We have furthermore investigated the relationship be- 
tween M imp at non-zero magnetic field and the impurity 
contribution to the susceptibility Ximp at zero field. At 
low temperature fc^Tc \D\, Ximp obeys a Curie law with 
a Curie constant interpolating between the free isotropic 
value of S(S + l)/3 for \D\ — > and the free anisotropic 
low temperature value of S 2 for \D\ — > 00. 56 It turns out 
that there is a simple relation between the Curie constant 
and the zero-temperature magnetization M imp for small 
magnetic fields gs^BB <C \D\: 



kBTx 



imp 



fM imp \ 



B=0, k B T<^\D\ 



W r ' h<\D\, fc B T<h 

(17) 

with h = gs^BB. The relative deviation between left 
hand and right hand side of Eq. (17), as determined by 



NRG calculations for all parameter combinations used 
in Fig. [4j is less than 1 %o. The relationship between 
zero-field susceptibility and magnetization expressed by 
Eq. (JTrj) is actually the same as for a doublet composed 
of states with magnetic quantum numbers ±M. In par- 
ticular, a free spin with easy axis anisotropy effectively 
reduces to such a doublet at low temperature UbT <C \D\, 
as discussed at the beginning of this section. 

In summary, the following picture of the low- 
temperature properties of a Kondo impurity with easy 
axis anisotropy is obtained: The impurity effectively 
acts as a two-level system with residual entropy Simp = 
ks In 2,^3 that splits up in a small magnetic field (leading 
to Simp = 0) and displays an unusual groundstate magne- 
tization indicative of a so-called "fractional spin"^. De- 
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FIG. 6. (Color online) Energy eigenvalues with magnetic 
quantum numbers M (upper panels) and magnetization 
(lower panels) as function of magnetic field for an anisotropic 
spin, described by Hamiltonian |5]l with D > (hard axis 
anisotropy), with a) S = 1, b) S = 3/2, and c) S = 2. For 
the magnetization curves and fields larger than the respective 
saturation field, temperature increases from top to bottom. 

spite this special property, the combination of Ximp and 
Mi mp shows that the magnetic response at low tempera- 
ture and field still resembles that of an ordinary magnetic 
doublet. 



VI. IMPURITIES WITH HARD AXIS 
ANISOTROPY 

We now investigate how an additional hard axis 
anisotropy (D > 0) affects the magnetic properties of 
the impurity. To lay the foundations for a study of the 
full impurity problem, we first discuss the magnetic field 
dependence of the magnetization for a free anisotropic 
spin described by Hamiltonian ^ with D > 0. 

For positive D and B = 0, the eigenvalues of Hamilto- 
nian ([5]) are energetically ordered according to the abso- 
lute value of their respective magnetic quantum number 
M. Depending on the spin S, the groundstate is thus 
either a singlet with M = (for integer S) or a dou- 
blet with M = ±1/2 (for half-integer S). In either case, 
the rest of the energy spectrum is comprised of doublets 
with magnetic quantum numbers ±M and 1/2 < M < S. 
The energy gap A| M | between a level with quantum num- 
ber M and the next higher-lying doublet is given by 
A| M | = (2\M I + 1)D. As a consequence of the magnetic 
field dependence of the eigenvalues described by the Zee- 
man term in Eq. |5]), n groundstate level crossings oc- 
cur for positive magnetic fields, with n = S for integer 
spin and n = S — 1/2 for half- integer spin. At the field 
Bm = &-\m\I Qs^b, the magnetic quantum number of the 
groundstate abruptly changes from — M to — (M+l) and 
hence the zero-temperature magnetization curve displays 
a discontinuous step. This effect is illustrated in Fig. [6] 
for spin S = 1,3/2, and 2. Finite temperature smears 
out the magnetization steps and renders them continu- 



ous. As the low-energy situation is the same in the vicin- 
ity of each groundstate level crossing, so is the effect of 
moderate temperature (cf. Fig. |6]c). 



A. Magnetic field dependence of the impurity 
magnetization 

We begin with the discussion of the magnetic field de- 
pendence of the impurity magnetization Ai{B) for equal 
g-factors and quasi-vanishing temperature T « 0. Mag- 
netization curves for impurity spin S = 1,3/2, and 2, 
and several values of the anisotropy parameter D > 
are shown in Fig. [7] Since the coupling strength pj , and 
thus the isotropic energy scale ksTn according to Table 
|TJ is kept constant, the ratio ksTn/D is varied. Note 
that a linear magnetic field scale is used in Fig. [7] to 
allow for an easy comparison with the results for a free 
anisotropic spin from Fig. [6] 

For an interpretation of the results for Ai(B,T ss 0), 
let us first consider the two limiting cases in which D 
is either small or large compared to UbTh- In the fol- 
lowing, imagine that we move from large magnetic fields 
to lower fields. If D is large, then little isotropic Kondo 
screening can occur before the anisotropy becomes ef- 
fective. As a guideline, we might thus think of the en- 
ergy spectrum of a free anisotropic spin S. On the other 
hand, for small D significant isotropic Kondo screening 
can take place before the magnetic field reaches the en- 
ergy scale defined by the anisotropy, so that it eventu- 
ally becomes more appropriate to think of the energy 
levels of a free anisotropic spin S — 1/2. As illustrated 
in Fig. |6j groundstate level crossings occur for a spin 
with hard axis anisotropy at certain fields and give rise 
to steps in the zero-temperature magnetization curve. 
For small and large anisotropy D, the impurity mag- 
netization J\4(B,T ss 0) also displays sharp, step-like 
features that are surrounded by magnetic field domains 
in which M. only slowly increases with B ("pseudo- 
plateaus"). However, due to the energy continuum of 
electronic states, these sharp features stay continuous 
even in the limit of zero temperature. The number of 
steps and their position relative to the field Djgs^B de- 
pends on the impurity spin S for large D and on the 
residual spin S ~ 1/2 for small D, respectively. The case 
S = 1 is special because the residual spin-1/2 cannot 
have uniaxial anisotropy. The single step which exists 
for S = 1 at large D therefore disappears for smaller val- 
ues of D. With respect to the energy scale imposed by 
the anisotropy, the steps in the impurity magnetization 
become well-defined for small and large D. In addition, 
the pseudo-plateaus become flatter. Fig. [7]c furthermore 
suggests that the two steps appearing in M(B,T w 0) 
for S = 2 and large D have different width. We are going 
to discuss the aspect of step width in more detail in Sec. 
|VI C| In particular, it will be shown that an impurity 
magnetization step is steeper if it occurs at larger field. 
It turns out that a standard z-averaging of the NRG re- 
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FIG. 7. (Color online) Impurity magnetization A4 for differ- 
ent anisotropy parameters D > as function of magnetic field 
for temperature ksT/W « 1.54- 10 _lj ~ 0, coupling strength 
pj = 0.07, equal g-factors, and impurity spin a) S = 1, b) 
S = 3/2, and c) S = 2. The value of D increases from bottom 
to top or from left to right, respectively. Note the use of a 
linear magnetic field scale and the rescaling of the Zeeman 
energy gsPBB by D. 



suits introduces artifacts into the magnetization steps for 
large anisotropy D. This problem is investigated in more 
detail in the context of an effective model at the end of 
Sec. |VI C| The plots shown in Fig. [7] are not visibly 
affected by this numerical shortcoming. 

The position of the steps in the impurity magnetiza- 



tion curves for both small and large anisotropy D seems 
interesting. Figs. [JJb and c show that for small D and 
impurity spin 5 — 3/2 and 5 = 2 a step occurs at a 
magnetic field which is smaller than the corresponding 
level crossing field for a free anisotropic spin 5—1/2. In 
contrast, for large D and all three impurity spins con- 
sidered in Fig. [7j each impurity magnetization step is 
found at a field exceeding the corresponding level cross- 
ing field for a spin 5. One might wonder whether the 
half-bandwidth W of the electrons has an impact on these 
two effects. To investigate this question for the case of 
small anisotropy, we have calculated additional magneti- 
zation curves for impurity spin 5 = 3/2 and 5 = 2 with 
decreasing coupling strength pj (0.09, 0.07, and 0.05). 
Since p = 1/2 W, a reduction of pj can be interpreted 
as an increase of the half-bandwidth W for fixed J. For 
each value of the coupling strength the anisotropy pa- 
rameter was chosen according to Table [T] so as to give 
a constant ratio UbTh/D = 1000. With this choice of 
D, a reduction of pj leads to a shift of the step in the 
impurity magnetization curve towards smaller fields rel- 
ative to D/gsp,g. This suggests that the effect observed 
for small D is not bandwidth-related. The question at 
which fields impurity magnetization steps occur for large 



anisotropy is considered in chapter VI C 



The limiting cases of small and large anisotropy D 
are connected by a regime with partial isotropic Kondo 
screening in which it is not possible to exclusively think 
in terms of the impurity spin 5 or an residual spin 
5 — 1/2. Impurity magnetization steps are broadened in 
this regime with respect to the energy scale D and, upon 
reducing the anisotropy parameter, move towards lower 
fields relative to D/gsptB (see Fig. [jjb). In addition, one 
step disappears for D ss ksTn in case of integer impurity 
spin 5 (cf. Fig. [7] c) . Isotropic Kondo screening effec- 
tively changes the impurity spin 5 from integer to half- 
integer and vice versaP*l Depending on the ratio UbTh / D 
and the resulting degree of isotropic Kondo screening, 
different behavior of the impurity magnetization is there- 
fore observed, in particular for fields gsPBB < D: While 
there is little magnetic response for an effective integer 
spin (at zero temperature, there is none at all for a free 
anisotropic integer spin as shown in Fig. [6]), we have 
larger impurity magnetization M. ~ for effective 

half-integer spin. However, in the latter case an addi- 
tional pseudo-spin- 1/2 Kondo effect occurs^ which leads 
to a suppression of the impurity magnetization for mag- 
netic fields gspBB <C D. 

In order to study the three different anisotropy regimes 
{k B T H /D < 1, k B T H /D sa 1, and k B T H /D > 1) in 
more detail, selected impurity magnetization curves from 
Fig. [7] are replotted using a logarithmic magnetic field 
scale (see Figs. [8] [9j and 10). This brings out effec- 
tive spin- 1/2 Kondo effects more clearly and allows for a 
better comparison with the behavior of an isotropic im- 
purity. Figs. [8] [9j and [10] furthermore demonstrate the 
effect of finite temperature on the impurity magnetiza- 
tion for the case of hard axis anisotropy. 
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FIG. 8. (Color online) Impurity magnetization M as func- 
tion of magnetic field for impurity spin 5 = 1 and anisotropy 
parameter a) D/W = 1(T 2 , b) D/W = 1(T 6 , and c) D/W = 
10~ 7 (solid lines, cf. Fig. [7]). Dashed lines show M(B) for the 
isotropic case and k B T/W « 1.54 • 10" 15 « 0. For the light 
gray lines, the approximate value of ksT/W increases from 
left to right from a) 4.35 • 1(T 4 to 6.10 • 10~ 2 , b) 1.99 ■ 10~ 7 
to 6.79 ■ 1(T 3 , and c) 4.73 • 1(T 10 to 6.79 ■ 1CT 3 . Thin vertical 
lines indicate the respective Zeeman energy gs^BB — D. For 
the chosen coupling strength pj = 0.07 and D = 0, we have 
k B T H /W w 3.39 • 10" 7 according to Table |l| 
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FIG. 9. (Color online) Impurity magnetization M as function 
of magnetic field for impurity spin S — 3/2 and anisotropy 
parameter a) D/W = 10" 2 , b) D/W = 10~ 6 , and c) D/W = 
10~ 10 (solid lines, cf. Fig. As before, dashed lines show 
A^B) for the isotropic case and ksT/W « 0. For the light 
gray lines, the approximate value of ksT/W increases from 
left to right from a) 1.61- 10~ 5 to 1.05- 10"\ b) 3.45- 10~ 7 to 
6.79 • 10~ 3 , and c) 1.01 ■ 10" 11 to 6.79 • 10" 3 . Magnetic fields 
satisfying g$PBB — nD with n = 1,2 are highlighted by thin 
vertical lines. According to Table [j] k B T H /W w 4.55 ■ 10~ 7 
for pj = 0.07 and D = 0. 



At large magnetic fields gsHsB D, the magneti- 
zation curve A4(B,T w 0) for an anisotropic impurity 
closely resembles the isotropic result. For lower fields 
9s^bB > D, the anisotropy eventually becomes effec- 
tive and the impurity magnetization begins to deviate 
from the isotropic curve. In particular, A4(B,T s» 0) 
displays a linear dependence on B for magnetic fields 
that are small compared to all relevant energy scales. If 
temperature is high so that kgT ^> D, then Ai(B) is 
suppressed for magnetic fields of the order of and smaller 



than ksT / gsp^B- On the other hand, judging by the rel- 
ative deviation from the zero-temperature curve, finite 
temperature has negligible effect on the impurity magne- 
tization if the thermal energy falls into the energy regime 
in which A4(B, T w 0) oc B (such a temperature indepen- 
dence is known for an isotropic impurity with S = 1/2, 
as discussed in Sec. IV A). In this regard, M differs from 



the magnetization of a free anisotropic spin for which 
non-zero temperature is always relevant. If the impurity 
magnetization curve features steps, they are smeared out 
for sufficiently high temperature. In case there is more 
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FIG. 10. (Color online) Impurity magnetization M as func- 
tion of magnetic field for impurity spin S = 2 and anisotropy 
parameter a) D/W = 10~ 2 , b) D/W = 10~ 6 , and c) 
D/W = 10 -8 (solid lines, cf. Fig. u\. Dashed lines again 
show A4(B) for the isotropic case and ksT/W « 0. For the 
light gray lines, the approximate value of ksT/W increases 
from left to right from a) 7.54- 10" 4 to 1.05- 10 _1 , b) 3.45- 1CT 7 
to 6.79 • 10" 3 , and c) 5.26 ■ 10" 11 to 6.79 • 10" 3 . Thin ver- 
tical lines indicate magnetic fields satisfying gsfiBB — nD, 
with n = 1,2,3. For p.J — 0.07 and D = 0, we obtain 
k B T H /W w 6.8 • KT 7 . 



than one step in M.(B,T ~ 0) (cf. Fig. 10 a), then, as 
a further difference compared to a free anisotropic spin, 
non-zero temperature has unequal effect on the different 
steps. We come back to the last two observations in Sec. 

rvrbi 

Ref. [5(3 arrives at the conclusion that an impurity with 
additional hard axis anisotropy undergoes anisotropic 
spin-1/2 Kondo screening if higher-lying impurity states 
are frozen out so that the impurity is effectively reduced 
to a spin-1/2 doublet. For gsfiBB <C D and low tem- 



perature, the impurity magnetization is therefore either 
suppressed if the impurity is effectively reduced to a sin- 
glet (cf. Figs. [8] a, [9] c, and 10 a), or it reflects the 



magnetic response of a Kondo screened spin-1/2 doublet 
(cf. Figs. [8] c, [9] a, and [T()| c) . For small D and inte- 
ger impurity spin, the corresponding Kondo temperature 
Tft decreases faster than linear in D upon reducing the 
anisotropy (it even drops exponentially fast in the spe- 
cial case S = 1 j."' 1 ' As a result, the pseudo-plateau that 
appears for S — 2 at magnetic fields gsP^B < 2D (see 
Figs. 10 c and[7]c) becomes flatter and broader when de- 
creasing D, whereas a reduction of D eventually leads to 
quasi-isotropic behavior over the whole considered mag- 
netic field range for S = 1 (see Figs. [8]c and [7] a). On 
the other hand, for large and increasing anisotropy D and 
half-integer impurity spin, the observed screening effect 
is increasingly well described by an exchange-anisotropic 
S = 1/2-Kondo model, whose Kondo temperature has a 
value much smaller than W/ ks This means that the 
pseudo-plateau occuring for impurity spin 5 = 3/2 and 
fields gsHBB < 2D (see Figs. [9]a and[7]b) becomes more 
pronounced for larger anisotropy D. As a final remark, 
Fig. [TO] a once again shows the different width of the two 
impurity magnetization steps for S = 2 and large D. 



B. Influence of the electron g-factor on A4 and the 
connection between Ml and Mi mp 

We now investigate how the impurity magnetization 
M(B,T w 0) is affected by a non-zero electron g-factor 
corresponding to a positive ratio g e / 9s > 0- Because of 
the sharp features that are found in the magnetization 
curves, the case of hard axis anisotropy seems well suited 
to study the influence of non-zero g e (alternatively one 
could examine the effect of the electron g-factor on the 
basis of the linear magnetic field dependence of Mi(B) for 
small B). As an example, Fig. 11 shows magnetization 
curves Ai{B, T « 0) for impurity spin S = 1, moderately 
large D, and several g-factor ratios interpolating between 
a local field (g e = 0) and a bulk field (g e = g$)- 

The results presented in Fig. [TT| demonstrate that a 
positive electron g-factor effectively causes a rescaling of 
the magnetic field: A ratio g e /gs > shifts the impurity 
magnetization curve for g e — towards larger fields and 
thus reduces A4(B,T w 0), which is a monotonically 
increasing function of B, for a fixed magnetic field value. 
Using the notation M.(B, T, gs, g e ), this statement can 
be expressed as: 



M(B,T^0,g s ,0) =M(B,T^0,g s ,g e ) , (18) 

9e 



B 



9s, 



(19) 



with a rescaling factor r](pj, g e /gs) ^ 1 for g e /gs ^ 
that depends on the coupling strength pj . Taking the 
magnetization curve for a local field as reference, we 
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FIG. 11. (Color online) Impurity magnetization M as func- 
tion of magnetic field for impurity spin S = 1, D/W = 10~ 4 , 
coupling strength pj — 0.07, and several values of the g-factor 
ratio g e /gs (cf. Fig. [7]a for g e — gs). For equal g-factors, the 
impurity contribution to the magnetization M; mp is shown, 
too. The temperature is k B T/W « 1.54 • 10" 15 « in all 
cases. 



lacks sharp features in the isotropic and easy axis case, 
the effect of a non-zero electron g-factor is more subtle, 
though. 

The numerical results thus strongly suggest that the 
rescaling factor r)(p.J, 1) and the proportionality factor 
a(pJ) take the same value. Furthermore, we find our 
results to be compatible with the conclusion of Ref. 05] 
that, in case of the isotropic Kondo model with 5 = 
1/2, the g-factor of the electrons (or equivalently their 
magnetic moment) is irrelevant for impurity properties in 
the limit of infinite bandwidth, corresponding to p J — > 0. 

In the case of impurity spin 5=1/2 and vanish- 
ing anisotropy D, it is also possible to compare the ob- 
tained values for r/(pj, 1) and a(pJ) with previously pub- 
lished results. For the exchange- anisotropic multichannel 
5 = 1/2-Kondo model with transverse coupling strength 
pj± <C 1 (see Eq. ( 29 1 for the meaning of the symbols 



J± and Jii), it is known that the impurity contributions 
to the free energy Fi mp (B,T,gs,g e ) for local and bulk 
magnetic field have the following relation!^ 



may therefore state that the impurity effectively "feels" 
a smaller magnetic field if there is also a Zeeman term 
for the electrons in Hamiltonian ([I]). 

For a bulk field, Fig. \TT\ additionally shows the impu- 
rity contribution to the magnetization M- lmp . Note that 
M and M imp are equal by definition for a local field . 
As previously found in the isotropic (cf. 
and easy axis case (cf. 



Sec. IV A) 



Sec. 



V A ) , the numerical re- 
sults strongly suggest a proportionality A4(B,T w 0) = 
a(pJ)M imp (B, T«0) for fields gsHnB < W. A study of 
the connection between M. and M- lrap for anisotropy pa- 
rameters < D -C W and impurity spin 5 = 1,3/2, and 
2 reveals that the average values for a(0.07) are nearly 
identical to the values for the isotropic case reported in 
Table |TJ and that the standard deviations have the same 
order of magnitude P3 In combination with the previous 
results for D < 0, we hence conclude that the effect of 
the uniaxial anisotropy on the value of a must be very 
small as long as \D\ -C W. 

An analysis of the magnetization curves for hard axis 
anisotropy furthermore shows that the values of the 
rescaling factor 77(0.07, 1) (approximately 1.038 for 5 = 1, 
1.0375 for 5 = 3/2, and 1.037 for 5 = 2) and the corre- 
sponding proportionality factor a(0.07) are remarkably 
similar. In case of 5 = 1, we have also studied the dif- 
ference between a local and bulk magnetic field for the 
two other coupling strengths previously considered (i.e., 
for pj = 0.05,0.09). The rescaling factors obtained for 
D < W (77(0.05,1) « 1.0265 and »?(0.09, 1) w 1.0497) 
are again in remarkable agreement with the correspond- 
ing values of a(pJ). Moreover, we find that Eqs. (18) 



and (19), with the values of rj(p,J,l) as determined for 
D > 0, are also suitable to describe the relation between 
the impurity magnetization curves for a local and bulk 
magnetic field for anisotropy D < 0. As Ai(B,T m 0) 



F imp (B,T, g Sl gs) 
B 



F imp {B,T,g s ,0) , (20) 
(l-2fS/n)B . (21) 



Here, / is the number of electron channels and 5 is 
the phase shift generated by the longitudinal coupling 
Jii . This result for i^imp is a generalization of the 
conclusion that the impurity contribution to the sus- 
ceptibility Ximp(dS)9e) a t zero magnetic field satisfies 
Xim P (.gs,3e) = Hg e /gs,S)xim P (gs,9s), with a certain 
factor A, for the single-channel exchange-anisotropic 5 = 
1/2-Kondo model with pj± -C 1.^ Using the defini- 
tion for the impurity contribution to the magnetization, 
-Mimp = —dFi mp /dB, and the equivalence of M; mp and 
the impurity magnetization A4 for g e — 0, the following 
relation is obtained from Eqs. (20) and (21): 



M imp (B,T,g s ,g s ) = 

(1 - 2/5/tt) M((l - 2f5/*)B, T, g s , 0) . (22) 

On the other hand, the proportionality A4(B,T pa 
0,9s, 9s) = a(pJ)M inip (B,T rj 0,gs,9s) implied by our 
NRG results can be combined with Eqs. ( 18 ) and ( 19 ) 
to give: 



M imp (B,Tn0,g s ,gs) 



1 m(- b 



,T w 0, 5 5,0 . (23) 



With / = 1 and the phase shift for the case of an elec- 
tron band of width 2W with constant DOS p = 1/2W, 76 
S(pJn) = arctan (7rpJ||/4) (note the sign change with re- 
spect to Ref. [7S]), we compare Eqs. (22) and p3| and 
deduce for pj± <C 1: 
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FIG. 12. (Color online) Proportionality factor a in A4 = 
aMimp as function of the transverse coupling strength pj± 
for impurity spin S = 1/2, vanishing uniaxial anisotropy D, 
and three values of the longitudinal coupling pju . Dashed 
horizontal lines mark the values of a, given by the numbers 
in parentheses, that are predicted by Eq. ( |24| ). Open symbols 
indicate the proportionality factors for the isotropic case (i.e., 
J|| = Jx), which are also found in Table [I] As before, a has 
been averaged over magnetic fields 

for k B T/W « 1.54 ■ 1CT 15 w 0. The corresponding standard 
deviations would amount to error bars smaller than the sym- 
bol size. 



«Hi) = v(pJ\\A) 



1 - ^ arctan (7179 Jy/4) 



(24) 



This equation predicts, in particular, that both the pro- 
portionality factor and the rescaling factor tend to 1 in 
the limit pJ\\ — > 0. From Eq. (l24| the following values for 
a(pJ\\) are obtained: a = 1.02563 (pJn = 0.05), 1.03623 
(0.07), and 1.04704 (0.09). In Fig. pj we present NRG 
results for a for an impurity spin S — 1/2 with van- 
ishing uniaxial anisotropy, but with additional exchange 
anisotropy. It is seen that the calculated proportional- 
ity factors indeed approach the predictions of Eq. ( 24 ) 
for decreasing transverse coupling strength pJ±. For 
pj± = 0.01 and all three considered values of pj\\ , the 
relative deviation is about 4 • 10 -5 . In the isotropic case 
(i.e., J|| = Jj_ = J), the relative deviation is less than half 
a percent, with better agreement for smaller coupling J. 

Since Eqs. ([20| and (21) also hold for non-zero tem- 



perature, the connection between A4(B,T > 0) and 
M- imp (B,T > 0) has been studied for isotropic impu- 
rity spin S = 1/2,1, and 3/2, and isotropic coupling 
pj = 0.05, 0.07, and 0.09. It turns out that M and M imp 
are still proportional for non-zero temperature and that, 
for gspsB/W E [10~ 13 , 10 _1 ], the relative deviation be- 
tween the proportionality factor and the corresponding 
value a(pJ,T 0) is less than 1 %o for thermal energies 
k B T/W < 10~ 2 . 



C. Field-induced Kondo effect 

To better understand the steps in the low tempera- 
ture magnetization curves for large anisotropy D > 
(cf. Figs. [7] and [13 a), we now derive effective (simpli- 
fied) models near groundstate level crossings (LCs) of the 
corresponding free anisotropic spin (see Fig. [6]). These 
models are approximations to the full Hamiltonian in the 
limit D — > oo. 



1. Effective models near groundstate level crossings in the 
limit of arbitrarily large anisotropy 

For given impurity spin S > 1 let us consider one of 
the groundstate level crossings of the corresponding free 
anisotropic spin (cf. Fig. [6]). If D is large, then the 
two levels which cross in the groundstate are energeti- 
cally well separated from the rest of the spectrum in the 
vicinity of the LC field. As an approximation for the 
full impurity model near this free LC field, we therefore 
project the impurity degrees of freedom in Hamiltonian 
([I]) onto the two impurity states involved in the free LC. 
This way, the impurity Hilbert space is reduced to two 
states and the impurity spin S can thus be mapped to 
an effective spin-1 /2. While the projection becomes exact 
only in the limit D — ¥ oo, we expect it to be a quantita- 
tive approximation for D ^> W and D ^> fc#T. 

The mapping of the impurity spin to a spin-1/2 is an 
extension of the ideas from Rcfs. l6Tlandl56l In contrast 
to the case of zero magnetic field that has been stud- 
ied there, we do not project onto impurity doublets with 
M = ±1/2 (see below). Furthermore, at each LC, i.e., 
for each step in the magnetization curve, the impurity is 
reduced to a different pair of states. As a consequence, 
different parameters of the effective model are obtained 
at each LC. 

We intend to use the effective models to determine the 
magnetic fields at which steps appear in the impurity 
magnetization curves for large anisotropy D (cf. Fig. [7|), 
and to investigate how the properties of the full impurity 
model differ near the various free LCs (as indicated, e.g. 
by Fig. 



13 a). Compared to the full model, the effective 



models are numerically less demanding as they feature an 
impurity spin-1/2 independent of the value of 5 1 , and they 
allow to study the effect of the different terms appearing 
in the effective Hamiltonian. 

To be specific, we consider the two impurity states 
with magnetic quantum numbers — M and —{M + 1) 
(assuming M > 0), which cross at the free LC field 
Bm — (2M + 1)D / gsP-B, and project Hamiltonian 
onto them. The effective model is determined by requir- 
ing that its matrix representation be equal to that of the 
full model in the chosen subspace. Note that we have to 
introduce new impurity states by shifting the magnetic 
quantum number in order to map the impurity spin to a 
spin-1/2. This mapping then corresponds to the follow- 
ing replacements for the impurity spin operators: 
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S x -> y/(S-M)(S + M + l) s x , (25) 
S y -> y/{S-M)(S + M + l) s y , (26) 
S z -> s z -(M + l/2) l a . (27) 

The impurity now has spin s = 1/2 and the replacement 
for (S z ) 2 directly follows from that for S z . In the param- 
eter regime in which the projection is valid, the mapping 
for S z leads to the following connection between the im- 
purity magnetization of the full and effective model: 



M/g sf i B = ~(S Z ) « -(s>) + (M + 1/2) 



(28) 



There is an analogous relationship for the impurity con- 
tribution to the magnetization Mj mp . According to its 
definition, the impurity contribution to the magnetic sus- 
ceptibility Ximp is not affected by the shift of the magnetic 
quantum numbers. 



Applying mappings ( 25 ) to ( 27 ) to the full impurity 
model and dropping all constant terms, the following 
Hamiltonian is obtained: 



H S (S, M)=J2 £ kacl a C ka - KS Z 



(29) 



+ 9s^b{B - B M )s z , 



with the set of parameters 



Efe 



Efc + Vge^BB 

(a = ±1/2) , 
K = (M + 1/2) J , 



J ± = y/(S — M)(S + M + 1) J 
J\\= J , 
B M = (2M + l)D/gsli B ■ 



(30) 

(31) 
(32) 
(33) 
(34) 



In contrast to the full Hamiltonian, H S (S, M) is exchange 
anisotropic with J± > J»: J±/J\\ grows with S and de- 
creases with increasing M (while always present, the ex- 
change anisotropy thus becomes weaker with every fur- 
ther LC). The Zeeman term for the impurity is now 
expressed relative to the free LC field Bm ■ H s (S, M) 
furthermore contains the new term — kSq representing 
an effective magnetic field, which couples to the elec- 
tron spin at the origin and points in the opposite direc- 
tion of the external field B. With respect to NRG, this 
term can be regarded as spin-dependent scattering at the 
zeroth site of the Wilson chain. It breaks the invari- 
ance under a spinflip transformation (cfc CT — > c^-a an d 



, — s z )), which H S (S,M) would otherwise 



possess for B = Bm. While the scattering parameter 



K grows with M, the ratio k/Bm is independent of M. 
Starting with the second LC, k is larger than Jy. The 
ratio k/J±, which at first is smaller than 1, also grows 
with M and eventually becomes greater than 1 if S is 
large enough. 

As an analogue to the free LC field Bm, we call the 
magnetic field -Belc f° r which the impurity magnetiza- 
tion vanishes at zero temperature (i.e., the two impurity 
levels are effectively degenerate), "effective level crossing 
(ELC) field": 



[s z )(B ELC ,T = 0) = 



(35) 



In the parameter regime in which the mapping to a spin- 
1/2 is valid, there is a step in the impurity magnetization 
curve of the full model at the ELC and, according to Eq. 
@, the value of M at the ELC is M « g s ^ B {M + 1/2). 



In the following, we discuss the properties of Hamiltonian 
( 29 1 in more detail for the two different cases g e > and 



9e = 0. 

Let us begin with the case g e > 0. As the free LC field 
Bm is proportional to D, the limit D — ¥ oo also corre- 
sponds to the limit B — > oo. A non-zero Zeeman coupling 
of the electrons therefore leads to their complete polariza- 
tion so that formally they may be replaced with spinless 
fermions (corresponding to spin-down electrons). Since 
the remaining fermion band is then completely filled, all 
interaction terms vanish and the electrons can be com- 
pletely eliminated from the problem. For g e > and 
arbitrarily large D, Hamiltonian (29) thus reduces to a 
pure spin model: 



>o) 



(B) = gs^B B 



gsVB J ~ 



(36) 



Here, we have introduced a relative magnetic field B = 
B — Bm- As the only remnant of the interaction be- 
tween impurity and electrons, a shift of the free LC field 
remains. This shift is positive for antiferromagnetic cou- 
pling J > and only depends on the coupling strength, 
but not on S or M. It is thus the same for all LCs. From 



the effective model (36) we learn that the ELC fields 



eventually exceed the free LC fields for g e > and large 
anisotropy D (cf. Fig. [7|). 

We now turn to the case of a local magnetic field. Set- 
ting g e = and using the relative field B, Hamiltonian 
(29) becomes the effective model for arbitrarily large D: 



H%-°\5;S,M) 



H S (S,M) 



9e=0 



(37) 



We 



(9e 



properties of 

: Bm,a — Bi\ 



are particularly interested in the 
eff =0) (B; S, M) at the ELC field B ELC = n ELC - n M - 
Due to the scattering term, the effective model does 
not exhibit a spinflip-invariance at the ELC. It there- 
fore seems that the ELC is not characterized by spe- 
cial symmetry properties. A spin- independent (poten- 
tial) scattering term can be treated by transforming to 
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scattering states which diagonalize the electronic part of 
the Hamiltonian (cf. App. C of Ref. [77j) . Although 
such a transformation can be easily adapted to the case 
of spin- dependent scattering, it does not seem to yield 
the intended results. The approximation which is used 
in the spin-independent case (a modification of the den- 
sity of states at the Fermi level expressed by an effective 
coupling parameter)-^ would restore spinflip-invariance 
for B = in the spin-dependent case and would thus 
erroneously imply (s z )(B = 0, T) = 0. Instead, we are 

going to use NRG to determine the ELC field -Belc an d 
to study the properties of ffgg e (-Belc! S, M). 

For the interpretation of the properties of the full 
Hamiltonian near the ELCs, we are going to use the 

main results for the effective model H^~ ^ (Belc! S, M). 
These are explicitly demonstrated in Sec. |VI C 4 and are 
summarized in the following. At an ELC the impurity 
spin s — 1/2 of the effective model is Kondo screened 
for T — > 0. The temperature dependence of the impu- 
rity contribution to the entropy at an ELC is described 
by the corresponding universal function for the isotropic 
S = 1/2-Kondo model with zero magnetic field. Since 
the parameters of the effective model arc different near 
each free LC (see Eqs. (31 ) and (32)), there is also a dif- 
ferent Kondo temperature Tk at each ELC. It turns out 
that Tk decreases with increasing M, i.e., Tk becomes 
smaller with every further ELC. 



2. Magnetic field dependence of impurity contributions near 
effective level crossing fields 

Armed with the effective model H^~ \B;S,M) for 
a local magnetic field, we now study in detail the field 
dependence of typical impurity contributions of the full 
impurity model for moderately large anisotropy D > 0. 
In Fig. [13] results for M imp (B), S imp (B), and T X im P {B) 
are shown for impurity spin S — 3 and anisotropy 
10 -3 . As before, equal g- factors have been 
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Let us start with a discussion of the magnetization 
curves depicted in Fig 
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FIG. 13. (Color online) Impurity contribution to a) the mag- 
netization Mi mp , b) the entropy Simp, and c) the effective 
moment Tximp as function of magnetic field for impurity spin 
S — 3, hard axis anisotropy D/W = 10 -3 , and several tem- 



13] a. According to the previ- perature values. In plot a, the dashed curve shows M imp for 



ously considered behavior of the impurity magnetization 
M. and its connection with Mj mp , there are also steps 

in M imp (B) at low temperature. These steps have finite (equal) peak positions for S imp and T X im P - As before, the 

0.07 and equal g-factors 



ksT/W ~ 0, and for the curves in plots b and c temperature 
increases from bottom to top. Thin vertical lines indicate the 



'imp 

widths for T — > and are smeared out for sufficiently high 
temperature. It is noticeable that the steps have differ- 
ent widths: those occurring at larger magnetic field are 
steeper. Fig. [13] a furthermore indicates that the influ- 
ence of non-zero temperature is different for the different 
steps. In contrast, for a free anisotropic spin the steps in 
the magnetization become discontinuous for T — > and 
the effect of non-zero (small) temperature is the same for 
all of them (see Fig. [6|. In the chosen representation of 
Fig. [l3]a, the pseudo-plateaus between the steps become 
flatter in the direction of increasing magnetic field and 
approach the true plateaus of the free anisotropic spin 



coupling strength is chosen as pJ 
are used. 



from below for growing anisotropy D. 



The behavior of M- lvxp {B) as shown in Fig. 13 



be understood by considering the magnetic field depen- 
dence of Mj mp for the isotropic S = 1/2-Kondo model 
with g e — gs- In this case, M imp (B,T = 0) is de- 
scribed by a universal function f(x) with the varia ble 
x = g s ti B B/k B T H (and T H cx T K , cf. Sec. |IVA[ )P 
f(x) is linear in x for x <C 1 and thus the slope of 
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Mimp(B) for small field s is higher if the Kondo tem- 
perature is smaller J22MI This relation is also expressed 
by Wilson's definition of the Kondo temperature!^ The 
zero-field susceptibility for T <C Tk is a Tff-dependent 
constant: Xim P (S = 0, T <C Tk) = const. /Tk- Com- 
bined with the prediction of the effective model with 
g e = for the Kondo temperatures at the different ELCs, 
this explains why the steps in the magnetization curve of 
the full model have different widths at zero temperature. 
In case of the S = 1/2-Kondo model, temperature has to 
reach the scale of Tk to become relevant for the zero-field 
susceptibility.^ For this reason, thermal broadening of a 
step in M- lmp (B) begins at lower temperature if the step 
occurs at a later ELC with smaller Tk- Furthermore, 
away from an ELC the magnetization reaches values of 
the order of the respective saturation value for smaller 
magnetic fields (relative to the ELC field) if the Kondo 
temperature at the ELC is lower. It subsequently en- 
ters the regime of very slow growth towards saturation, 
which shows up in Fig. [13] a in the form of a pseudo- 
plateau. This also explains why pseudo-plateaus between 
later ELCs with smaller Tk are flatter. 

We now continue with a discussion of the magnetic 
field dependence of the impurity contribution to the en- 
tropy and the effective moment. Results for iSi mp (j3) and 
Tximp(B) at low temperature ksT < D are shown in Fig. 
[j~3|b and c, respectively. In both cases we observe peaks 
of varying height and width whose positions coincide with 
those of the steps in Mi mp (-B). We find that the peaks be- 
come both higher and narrower with every further ELC. 
If the temperature is reduced, the peak heights decline 
and at the same time, if T is not too low, the peaks be- 
come sharper. It is noticeable that there is a temperature 
below which the width of the first peak in both Si mp and 
T^imp varies only little as a function of T. 

At zero temperature both Sj mp and Tximp vanish for 
all magnetic fields. In case of the entropy, the reason is 
that the magnetic field either leads to a non-degenerate 
groundstate or it creates an effective impurity doublet 
which is then Kondo screened. The effective moment, 
on the other hand, has to go to zero since the slope of 
M- imp (B) at zero temperature, i.e., Xi mp (-B,T = 0), is 
finite for all fields. For large anisotropy D, the tem- 
perature dependence of the peak heights is determined 
by the pseudo-spin-1/2 Kondo effects that take place 
at the ELCs (see the next section and Fig. 14 for de- 
tails). By recollecting results for the S = l72-Kondo 
modelpS we can furthermore understand the different 
peak widths and their temperature dependence. In case 
of the S — 1/2-Kondo effect, the Zeeman energy gsjiBB 
has to reach the energy scale of max (fcsT, &bTr-) in or- 
der to considerably suppress S lmv and Tximp^I 11 par- 
ticular, since the lowest temperatures considered in Fig. 
13 are smaller than Tk at the first ELC (cf. the in- 
dicated temperature range in Fig. 14 a), temperatures 
are reached for which the thermal broadening of the first 
peak is small. 
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FIG. 14. (Color online) Impurity contribution to a) the en- 
tropy Simp, b) the effective moment T^mp, and c) the mag- 
netic susceptibility Ximp as function of temperature for impu- 
rity spin 5 = 3 and hard axis anisotropy D/W = 10~ 3 . The 
chosen non-zero magnetic fields correspond to the peak posi- 
tions for Simp and Tximp according to Fig. |13| Also shown 
are NRG results for vanishing coupling strength (dashed lines) 
and for the isotropic case (dash-dotted lines). Solid vertical 
lines delimit the temperature interval considered in Fig. |13| 
whereas dashed vertical lines indicate the thermal energy for 
which fcflT = D. In the insets, data from the main plots is 
presented for a reduced temperature range. 



Temperature dependence of impurity contributions at 
effective level crossing fields 



Fig. 

TXimpi 



14 shows the temperature dependence of S[ mp , 
and Ximp at the three ELCs that occur for the 
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example considered in Fig. 13 In case of entropy and 



effective moment it is thus demonstrated how the peak 
heights in Figs. [13] b and c decrease when the tempera- 
ture is lowered. 

Let us begin with a discussion of the results for 5i mp (T) 

0, T) for the isotropic 



shown in Fig 
impurity with 3 



14 



5'i m p(-B 



3 interpolates between the limiting 
values of k B \n7 (for T -> oo) and fc B ln6 (for T = 0), 
according to the screening of half a magnetic moment. 
In contrast, for the anisotropic impurity in zero field all 
higher-lying impurity levels are frozen out on the tem- 
perature scale k B T ~ D so that <Si mp quickly drops to 
zeroP^ At high temperatures fc^T 3> D, on the other 
hand, the anisotropic impurity behaves more and more 
like an isotropic impurity. For the three ELC fields, 
we observe qualitatively different behavior at low tem- 
perature k B T < D compared to the case of zero field: 
Starting at S- lmp ~ k B hi2, the effective impurity dou- 
blet, which is formed due to the magnetic field, under- 
goes Kondo screening with a value of Tk that decreases 
with every further ELC. While there is fair agreement 
with the results for a decoupled impurity at k B T > D, 
low temperature k B T -c D reveals that the ELC fields 
are not equal to the free LC fields. For this reason, all 
states except the respective non-degenerate groundstate 
are ultimately frozen out for T — > 0. 

We continue with a discussion of the temperature de- 
pendence of the effective moment Tximp as illustrated 
in Fig. [14] b. For the isotropic impurity with 5 = 3, 
k B Txim P {B — 0,T)/(gsfi B ) 2 obeys a Curie law at both 
high and low temperature and goes from S(S + l)/3 = 4 
(for T -> oo) to (S - 1/2) (5 + l/2)/3 = 35/12 (for 
T = 0). For the anisotropic impurity in zero magnetic 
field, on the other hand, the increasing thermal reduc- 
tion to a non-magnetic groundstate for k B T < D leads 
to a vanishing effective moment at zero temperature.^ 
At very high temperature k B T < W (not shown), we find 
the expected agreement between Tximp(B = 0) for the 
isotropic and anisotropic case. The effective moment of 
the decoupled anisotropic impurity first goes to about 1/4 
at k B T < D (according to the susceptibility of a doublet 
with AM = 1), but then quickly drops to zero because 
of the non-degenerate groundstate with good magnetic 
quantum number. For the anisotropic impurity at the 
ELC fields, we again observe Kondo screening for r^0 
with different Tk, starting with an effective moment of 
k B T X im P /(gsVB) 2 « 1/4 at k B T < D. 

Finally, we take a look at the results for the impurity 
contribution to the susceptibility Ximp(T) that are pre- 
sented in Fig. 



14 



c. According to its definition, Ximp is 
just the slope of Mj mp and thus directly yields informa- 
tion about the width of the steps in the magnetization 
curve shown in Fig. [13] a. For an isotropic impurity with 
S > 1, the zero-field susceptibility at low temperature is 
described by a Curie law and thus diverges for T — » 0. 
In contrast, Ximp(B = 0,T) for the anisotropic impurity 
has a maximum at k B T D and vanishes for zero tem- 
perature. Since the ELCs lie close to the free LCs in this 



example, the susceptibility for the decoupled anisotropic 
impurity displays a maximum at a thermal energy of the 
order of the level splitting before falling off at low temper- 
ature. At the ELCs, Xi mp (T) for the anisotropic impurity 
saturates at a finite value for T — > that increases with 
decreasing Kondo temperature. 

Recently, it has been demonstrated that a field-induced 
Kondo effect also occurs for Hamiltonian with easy 
axis anisotropy D < 0, additional transverse anisotropy 
E, and a local magnetic field aligned along the x-axisP^ 



4- Properties of the effective model for vanishing electron 
g-factor 

We now return to the effective model for g e = given 
by Hamiltonian (37) in order to study its properties in 



greater detail. The ELC field £>elc and the Kondo tem- 
perature Tk oi the ELC field are determined as func- 
tion of the parameters J\\ , J±, and k of the effective model 
or, respectively, as function of the parameters J, S, and 
M of the full Hamiltonian (according to Eqs. (31) to 



(33)). As a prerequisite, we have to figure out how to 



reliably extract these quantities from the NRG results. 

To determine -BelCi t ne impurity magnetization in 
units of gsl^B for the effective model, —(s z )(B), is cal- 
culated for low temperature k B T <c W. In the vicinity 
of an ELC, i.e., near its root, the impurity magnetization 
depends linearly on the (relative) magnetic field B. The 
root, which corresponds to -Belc at T = 0, can therefore 
be determined by performing a linear fit to the numerical 
data. However, the following complication arises: The 
position of the root of (s z )(B) depends on the value of 
the twist parameter z and thereby on the discretization 
of the electron band. On the contrary, a physically mean- 
ingful result for the ELC field should display only a weak 
dependence on the numerical parameters of an NRG cal- 
culation in order to accurately reflect the continuum limit 
A — > 1. It turns out that a standard z-averaging, i.e., an 
averaging of the impurity magnetization curves for dif- 
ferent values of z at fixed temperature, is not reasonable 
at this point. Near an ELC, such an averaging in gen- 
eral introduces artifacts into the averaged curve because 
of non-linear components which some of the z-dependent 
curves might already comprise. Similar numerical arti- 
facts are found in the z-averaged magnetization curves 
of the full model for large hard axis anisotropy. Upon 
closer inspection, one discovers that z-averaging divides 
the total height of a magnetization step into smaller "sub- 
steps" of equal height whose number corresponds to the 
number of z-values used. 

For one set of parameters of the effective model, the de- 
pendence of the impurity magnetization root on the dis- 
cretization of the electron band is demonstrated in Fig. 
|15| There, the three common discretization schemes are 
compared for three different values of the discretization 
parameter A and, only in this example, for 16 values of 
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FIG. 15. (Color online) Impurity magnetization (in units of 
ffs^s) for the effective Ham ilto nian wit h g e = and pa- 
rameters according to Eqs. (311 to (331 as function of the 
relative magnetic field for ksT/W » 10~ 15 » 0. The re- 
sults have been calculated using the discretization schemes 
by a) Yoshida et al. (with correction factoi 113 1 14 1 79 1 A A ) W^1 
b) Campo and Oliveira} 211 and c) Zitko and PruschkeP 1131 
For each discretization scheme, results are presented for three 
values of the discretization parameter A and 16 values of the 
twist parameter z (i.e., Zi = i/16 with i G {1, 2, 16}). In 
plot c, data points are vertically offset to enhance legibility. 
Numbers in parentheses denote the ELC field Belc and the 
Kondo temperature T^ LC at the ELC field, respectively (cf. 
main text). 



the twist parameter z. In all cases, we observe a spread 
of the position of the impurity magnetization root with 
respect to z. This variation decreases for smaller values 
of A and is always largest when using the discretization 
scheme by Zitko and Pruschke (ZP). The spread due to z 
defines a magnetic field interval which, in case of the dis- 
cretization by Yoshida et al. (Y) and Campo and Oliveira 
(CO), moves towards larger fields when A is reduced. In 
contrast, the ZP discretization leads to nested intervals 
so that an interval for smaller A is wholly contained in an 
interval for larger A. In the special case z — 1 (which cor- 
responds to the smallest root), the CO and ZP discretiza- 
tion give the same resultP^ Note that the A-dependence 



of the data shown in Fig. 15 is consistent with an agree- 
ment of the results of all three discretization schemes in 
the continuum limit A — > 1. However, in order to obtain 
reliable information about the continuum limit it would 
be necessary to perform an impractical extrapolation in 
A when using the Y or CO discretization. On the con- 
trary, the ZP discretization apparently allows to make 
a dependable statement about the limit A — > 1 on the 
basis of results for only a single value of the discretiza- 
tion parameter: Fig. [15] suggests that the continuum 
value of -Belc nes m t ne magnetic field interval that is 
spanned by the z-dependent roots for A > 1. It turns 



out that, using the ZP discretization, one can obtain a 
better approximation for the ELC field by averaging over 
the z-dependent impurity magnetization roots since the 
resulting mean value displays only a weak dependence 
on A (cf. the numbers in parentheses in Fig. 151. The 



spread of the roots with respect to z then provides a safe 
error estimate for the mean value (amounting to a rela- 
tive deviation of about 3 to 4 % for A = 3). However, 
the dependence of the mean value on A indicates that 
such an error estimate is far too pessimistic. For the ZP 
discretization with A = 3 and four z-values, we expect 
that the relative error of the obtained ELC field £?elc is 
about one order of magnitude smaller than suggested by 
the roots' dependency on z. 

Having determined B^ LC for all values of z, we can 
study the thermodynamic properties of the effective 
model at the ELC field by calculating the impurity con- 
tribution to the entropy S? (T, B^ LC ) for each z at 
the respective ELC field. For thermal energies that are 
small compared to the bandwidth, the NRG results for 
Simp(T) -^elc) can be aligned with the known universal 
temperature dependence of the entropy for the isotropic 
Kondo model with S = 1/2. We find, however, that there 
is a dip in the entropy for thermal energies close to the 
band edge that becomes more pronounced for stronger 
scattering k. The Kondo temperature j^ c > z character- 
izing the universal behavior of Sf mp {T, -B EL c) is obtained 
in the following way for each value of z: By comparing 
with the Bethe Ansatz solution for the impurity contri- 
bution to the magnetization in Sec. |IV A[ the value of 
Tx (as occurring in the BA solution of the Kondo model) 
for the isotropic NRG results is known (cf. Table [I] and 
Eq. Il0|). As a first step, the result for 5?, p (T, B% hC ) is 
restricted to the linear low-temperature regime in which 
a continuous curve is produced using a linear fit. This 
fit then allows to determine the value of y^ LC ^ z by com- 
paring with an isotropic NRG result for S = 1/2 with 
known Kondo temperature T^°. A better approxima- 
tion for the Kondo temperature at the ELC field T^ LC is 
again obtained by averaging over the z-dependent values. 
Regarding the variation with respect to z and the error 
estimate for the mean value, comparable statements hold 
true as in the case of the ELC field. 

Before continuing, we note that the root of the im- 
purity magnetization depends on temperature. This is 
possible since its value is not determined by symmetry 
properties of the effective model. When the temperature 
exceeds the Kondo temperature at the ELC field, two 
effects eventually occur: The slope of the magnetization 
curve decreases and the root moves towards larger rela- 
tive magnetic fields. 

We now study how the parameter k, which corresponds 
to spin-dependent scattering at the zeroth site of the Wil- 
son chain with on-site energies eo^ = — eof > as seen 

from Eqs. (29) and (31), affects the values of -Belc an d 
T^ LC . To this end, we interpret free parameter 

of the effective model. For the isotropic Kondo model 
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FIG. 16. (Color online) Kondo temperature at the ELC field 
for the effective Hamiltonian with g e = as function of the 
scattering parameter k (corresponding to a pair of on-site 
energies £04. = — e ot > for the zeroth site of the Wilson 
chain) for fixed J11 and three values of the coupling parameter 
J± > J|| . Numbers in parentheses denote the respective ELC 
field and lines are intended as a guide to the eye. Pessimistic 
error bars would be smaller than the symbol size. Additional 
crosses indicate the Kondo temperature for parameters which, 
according to Eqs. (311 and (321, correspond to the given 



values of the quantum numbers S and M. Small dots mark Tk 
for the case of potential scattering, i.e., for spin-independent 
on-site energies at zero magnetic field. 



in zero magnetic field it is known that ordinary (spin- 
independent) potential scattering can be approximately 
mapped to a modified electron DOS at the Fermi level 
or, equivalently, to an effective coupling parameter J e ff=^ 
This approximation predicts that an increase of the scat- 
tering parameter reduces J c g and thus also the Kondo 
temperature. 

In the example shown in Fig. |16| spin quantum num- 
ber 5 = 3 (as in Fig. 13 1 is considered and those cou- 



pling parameters J± are chosen which, as per Eq. ( 32 1 , 



are assigned to the three magnetic quantum numbers al- 
lowed for this value of 5 (i.e., M = 0,1,2). Without 
scattering ter m, th e Kondo temperature decreases when 
J_i_ is reducedpESl Additional spin-dependent scattering 
further lowers T^ LC just as standard potential scatter- 
ing with eo^ = eot at zero magnetic field does, but in 
comparison leads to smaller values of the Kondo temper- 
ature. In accordance with the expression for J c g from 
Ref. 1771 the sign of the spin- independent on-site energies 
does not affect Tk- Fig. [16] reveals that the decrease of 
T k AjC accelerates with growing scattering strength. Fur- 
thermore, we observe that the spin-dependent scattering 
has a larger influence on the Kondo temperature at the 
ELC field when the coupling parameter J± is smaller. 

Let us now turn to the effect of k on the position of 
the ELC field. An additional spin-dependent scattering 
term breaks the spinflip-invariance and is therefore the 
very reason for a non-zero value of -Belo It thus seems 
plausible that a larger value of k also leads to a larger 
absolute value of the ELC field (cf. the numbers in paren- 



theses in Fig. 16). This increase of |Selc| decelerates 
with growing scattering strength. A closer look at the 
data reveals that K again has a stronger effect when the 
coupling parameter J± is smaller. 

Additional crosses in Fig. [16] mark the Kondo temper- 
ature for those values of the scattering parameter k that 
follow from Eq. ( 31 ) for the three considered M quantum 



numbers. We observe that the effective model predicts a 
decrease of T^ LC with growing M, i.e., with every fur- 
ther ELC. As the example demonstrates, this decline of 
the Kondo temperature is due to three cooperating ef- 
fects: 1) According to Eq. (32), J± becomes smaller 
when M is increased. 2) Simultaneously, k becomes 
larger. 3) Because of the decreasing value of J±, the scat- 
tering parameter additionally gains in importance. For 
the ELC fields gs^BBELc/W that belong to the three 
special values of T^ LC , we obtain the following results 
(the error estimates indicate the variance with respect to 
z): -6.59±g;^ • ltr 3 (S = 3, M = 0), — 1.78±g;g? • 10" " 



-0.25 

(M = 1), and -2.36l£o9 • 1CT 2 (M = 2). 

Finally, we investigate how the ELC field and the 
Kondo temperature at the ELC field depend on the pa- 
rameters of the full Hamiltonian, i.e., on J, 5, and M, 
with the parameters of the effective model given by Eqs. 
(31) to (33). First of all, we note that all obtained (rel- 
ative) ELC fields are negative. This supports the con- 
clusion that, in the impurity magnetization curves for 
equal g-factors and large hard axis anisotropy presented 
in Fig. [7] the free LCs are only exceeded because of 
the electrons' non-zero magnetic coupling. For 5 = 1, 
3/2, 2 and J/VF = 0.14 (as in Fig. [7}, the following 
values for gs/J-BB^hc/W are obtained: -3.90±g;^ • 10~ 3 
(5 = l,M = 0), -8.27^30 • 10~ 3 (5 = 3/2, M = 1/2), 



-4.95 



+0.14 



10- 



(5 = 2,M = 0), -1.31+^ 



4-0.04 
05 



10" 



-0.18 

(5 = 2,M = 1). 

In contrast to the prediction of the effective model with 
g e > from Eq. ( 36 1 , the ELC fields depend on the quan- 



tum numbers 5 and M for vanishing electron g-factor 
(see Fig. 17). With increasing value of 5, i.e., with in- 



creasing coupling parameter J±, the absolute value of 



the ELC field grows as already seen in Fig. 16 A larger 
coupling J increases all parameters of the effective model 
(■^11 > J±> and K ) and, as demonstrated by Fig. 17 thereby 

leads to a larger value of |-Belc|- It is furthermore evi- 
dent that, with growing J, the quantum number 5 gains 



in importance: Fig. 17 shows that the "curves" for fixed 
coupling strength "fan out" for larger value of J. The 
dependence of the ELC field on the magnetic quantum 
number is the result of two counteracting effects: 1) A 
larger value of M leads to smaller coupling strength J± 
which, on its own, would lower |-E>elc|- 2) On the other 
hand, the scattering term becomes stronger with increas- 
ing M and would, on its own, enlarge |Selc|- For the 
parameters considered in Fig. [17] there is a growth of 
|-Belc| with M for S < 7/2 that decelerates with in- 
creasing M. In the case of S = 4 and both J/W = 0.16 
and J/W = 0.18, we observe a decrease of the absolute 
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FIG. 17. (Color online) Negative value of the ELC field for 
a set of parameter combinations of the effective Hamiltonian 
with g e = which, according to Eqs. ( |31[ ) to ( |33[ ), correspond 
to certain spin quantum numbers S, magnetic quantum num- 
bers < M < S, and coupling parameters J. Projections 
of the data points onto two planes are shown as small dots 
and lines connect points with the same value of S and J. 
Pessimistic error bars would be smaller than the symbol size. 



value of the ELC field in the last step. 

To conclude this section, the Kondo temperatures be- 
longing to the ELC fields shown in Fig. [17] are pre- 
sented in Fig. [THJ As the main result we find that, 
according to the above explanation, the value of T^ LC 
increasingly drops with growing M. On the other hand, 
a larger value of S increases the coupling strength Jj_ 
For the parameters considered in 
the "last" value of the Kondo temperature (i.e., 
for M = S — 1 at fixed J) also grows with S. It turns 



and thus also Tf LC 



Fig. 18 



out that the influence of S on J^ LC is reduced for larger 
coupling J. As a consequence, the "curves" for fixed 
J are "focussed" in the direction of increasing coupling 
strength. We find that the dependence of the Kondo 
temperature on J is the result of two counteracting ef- 
fects: 1) Both coupling parameters Jn and J± grow with 
J and would, on their own, lead to a larger value of T^ LC . 
2) However, the scattering parameter k is also increased 
and would, on its own, lower the Kondo temperature. 
For all parameter combinations presented in Fig. [18 



Tj| LC is a monotonously increasing function of the cou- 
pling strength J. As a final observation, the relative 
decrease of the Kondo temperature between two consec- 
utive values of M (for fixed J) becomes smaller for larger 
quantum number S. 



5. Comparison of anisotropy- and field-induced 
pseudo-spin- 1/2 Kondo effect for half-integer impurity spin 

In the case of half- integer impurity spin S > 3/2, two 
different types of pseudo-spin- 1/2 Kondo effects occur for 
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FIG. 18. (Color online) Kondo temperature at the respective 
ELC field for the same set of parameter combinations of the 
effective Hamiltonian with g e = as in Fig. |17| Again, 
pessimistic error bars would be smaller than the symbol size. 
Three small crosses in the plane spanned by fcsT^ LC jW and 
J/W indicate the value of the Kondo temperature for the 
isotropic case without scattering term (as inferred from Table 



I). 



large hard axis anisotropy D. At zero magnetic field, the 
anisotropy splits up the impurity multiplet and a doublet 
with magnetic quantum numbers M = ±1/2 remains as 
the lowest-lying impurity level (cf. Fig. [6] b) . If the en- 
ergy gap to the impurity states with M = ±3/2, which 
is equal to 2D, is sufficiently large, this doublet under- 
goes spin- 1/2 Kondo screening (let us call this Kondo ef- 
fect "anisotropy- induced" ) At non-zero magnetic field, 
on the other hand, ELCs with associated "field-induced" 
Kondo effects occur, as discussed in the previous section. 
Both types of Kondo effects show up in the corresponding 
impurity magnetization curves. In Fig. [9] a, for example, 
we observe the magnetic response of the effective impu- 
rity doublet at low magnetic fields gs^sB/D <C 1, which 
is then followed by a step at gs^sB/D m 2 due to the 
ELC. 

One might wonder how the Kondo temperatures of 
the anisotropy- and field-induced Kondo effects com- 
pare. It has been shown in Ref. [SUthat the anisotropy- 
induced Kondo screening in t he limit D/W — ¥ oo is ex- 
plained by Hamiltonian (37) with B = and param- 
eters J|| = J, Jf = (S + 1/2)J, and k = 0. On 
the other hand, the pseudo-spin- 1/2 Kondo effect at the 
first ELC (i.e., M = 1/2) is described by the same 
Hamiltonian with B = -E>elc an d parameters J|| = J, 
jf LC = y/(S — l/2)(5 + 3/2) J, and k = J. In par- 
ticular, we thus have jf LC < Jf for S > 3/2. Since 
Fig. [16] shows that additional scattering k reduces the 
Kondo temperature, we can conclude that for the 
anisotropy-induced Kondo effect is always larger than 
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Tf LC at the first ELC and, according to Fig. [l8j at 
all following ELCs, too. 



VII. SUMMARY 

In this article, we have reported on Numerical Renor- 
malization Group (NRG) calculations for a Kondo model 
with additional uniaxial anisotropy D. Results have been 
presented for non-zero magnetic field B and different ra- 
tios g e /gs °f electron and impurity g-factor. 

For a bulk field (i.e., for equal g-factors), a comparison 
of low-temperature NRG results for the impurity mag- 
netization A4(B) and the impurity contribution to the 
magnetization Mi mp (B) reveals that M. = aAii mp . The 
proportionality factor a > 1 decreases for smaller cou- 
plings pjy and pJj_. Compared to the case of a local 
field (i.e., g e = 0), a non-zero electron g-factor effectively 
rescales the magnetic field argument of the impurity mag- 
netization AA(B) at low temperature. Calculations for 
isotropic exchange interaction (i.e., Jy = J± = J) sug- 
gest that both effects must have a weak dependence on 
D as long as the uniaxial anisotropy is small compared to 
the bandwidth. In addition, the corresponding values of 
the rescaling factor rj(g e / 9s = 1) and the proportionality 
factor a apparently coincide. 

For an isotropic impurity (D = 0), we find that M(T ss 
0, B), unlike M; mp , does not display universal behavior 
in the usual sense as already noticed by Lowenstein.^l 

With additional easy axis anisotropy (D < 0), the zero- 
temperature curve M- lrap (B) approaches a Independent 
value for small, but non-zero magnetic field, which cor- 
responds to the effective moment of the respective "frac- 
tional spin"^ as given by Txim P (B — 0). The magnetic 
response at small fields and low temperature (relative to 
1 1? |) is thus reminiscent of an ordinary magnetic doublet. 
Appropriately, the impurity magnetization M is well de- 
scribed by a rescaled and shifted Brillouin function in 
this regime. 

In the case of hard axis anisotropy (D > 0), a non- 
zero magnetic field can lead to "effective level crossings" 
(ELCs), at which pseudo-spin-1/2 Kondo screening oc- 
curs. For g e = and D/W — > oo, these field- induced 
Kondo effects are described by an exchange-anisotropic 
spin-1/2 Kondo model with additional spin-dependent 
scattering at the zeroth site of the Wilson chain. At 
the respective ELC field, this scattering leads to a reduc- 
tion of the Kondo temperature T^ LC in a similar way as 
ordinary potential scattering does at zero magnetic field. 
In particular, the effective model predicts that T^ LC de- 
creases with every further ELC. This agrees with the ob- 
servation that the steps in the magnetization curves for 
large D, which are due to the field-induced Kondo effects, 
become steeper in the direction of increasing magnetic 
field. 
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Appendix A: Numerical Renormalization Group 
calculations with conduction electron Zeeman term 

In this appendix, we briefly describe the changes to the 
standard NRG procedure^ which are necessary in order 
to carry out calculations with an additional Zeeman term 
for the conduction electrons. 



1. Logarithmic discretization 

The starting point is the continuous energy represen- 
tation of the Hamiltonian from Eq. ([7| with a restriction 
to the physically reasonable case h < W. In the following 
it is assumed that the magnetic field, which appears in 
h = g e p,BB, is non-zero and fixed and, to simplify the no- 
tation, that the chemical potential is zero. We introduce 
abbreviations for the absolute value of the integration 
boundaries in Eq. 



B± = \±W + ph\, 



(Al) 



rescale the integration variable e, and change to rescaled 
electron operators (cf. Ref. IT5]) : 



(A2) 
(A3) 

(A4) 
(A5) 




a ifi = \J Bp a efi for e < 



Using J_ w de p(e) = 1 and defining the normalized ze- 
roth state of the Wilson chain as 



J d£+ y/p(Z+B+-»h)B, 



(A6) 
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we then obtain an equivalent expression for the electronic 
and interaction term in Eq. Q: 



H, 



cb+int 



GsH&L ( A7 ) 



W 



Next, the logarithmic discretization of the conduction 
band is carried out according to one of the available dis- 
cretization schemes ^ 3 * 64 * 79 * ^ by dividing the integration 
range [—1,1] into standard intervals Zjjj and using the fol- 
lowing weight function on the m th positive and negative 
interval, respectively: 



^Pmfx (^m ) 



With a = ±, 



H Op 



y,h) 



ph) 



(A8) 



At this point, the substitution (All I is still valid for 



arbitrary pis). The above expressions simplify in the case 
of a constant density of states, p(e) = 1/2W, as a shifted 
constant DOS is, of course, still a constant DOS: 



B 



nifi 



H ch+int -^W ^ ^£^(h = 0) a s X a 



nifi 



(A12) 
(A13) 



Here, £ ^ and 7 m are the "energies" and expansion coeffi- 
cients, respectively, for the system with a local magnetic 
field (i.e., with g e =0). 



2. Tridiagonalization 



j ^ d^p(^-ph)W 



(A9) 



and new operators corresponding to the weight func- 



tions tPmfiiifi) on the intervals J^, we have the following 
exact expansion for the zeroth state of the Wilson chain: 



fan = % 



(A10) 



In addition, a dimensionless "energy" £^ has to be as- 
signed to each interval ZjJ, for each spin projection /i. This 
is done according to the chosen discretization scheme by 
using the weight function (A8) with the shifted DOS, 



leading to a discrete approximation to Hamiltonian ( A7 1 



H 



B s 

< l+int ->• W ^ 7U f mM gmp Sm^ 



(All) 



Since the rescaling factors 



B s 
n 

W 



P~ 



9e 9sPbB 



9s 



W 



(A14) 



depend on spin projection fi and magnetic field B, the 
tridiagonalization of Hamiltonian (A12), which leads to 
the Wilson chain with hopping parameters U^B) and on- 
site energies e^B), has to be done separately for spin- 
up and spin-down and for each value of B. In case of 
a constant DOS, Eqs. 1A.121 and (IA13I) show that the 



only necessary modification of an existing code solving 
the recursion relations given in Ref. [2] is to multiply all 
"energi es" £ ^ and coefficients 7^ with the appropriate 
factor jAlil. 



For a particle-hole symmetric DOS we have = 



-£ m -n and 7 

{-l) n v. 



infi 

and v„ 



In 



Using the Ansatz 



(-1)' 



infi 



for the coefficients 



of the orthogonal transformation (following the notation 
of Ref. IH]) . it can then be shown that t^(B) = Ui(B) 
and e^(B) 



-€ii(B) for all sites i of the Wilson chain. 
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